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Abstract 


A deeper understanding in the transverse impact of wire ropes is 
required to be able to design more and improved structures to reduce the 
soaring loss of human lifes in traffic accident, by fall of object in mines and 
the huge amount of rockfalls and landslides ubiquitous all along the rugged 
geography that we have. By all these reasons, in this thesis are given the 
mathematical models of the forces in the ends of a steel wire rope and the 
movement of the body that impact transversely in its middle point. To achieve 
these objectives, both ends of the wire rope are considered rigidly fixed, 
reducing its geometry to an unidimensional elastic string and modeling the 
body as a material point that could be a rock, automobile or any machinery 
that could fall over a protection barrier, achieving this way a high level of 
acceleration what in turn induces large displacements and internal stresses 
around the contact zone in the wire rope. Next the models are validated 
using the developed equations to calculate the reaction of the steel wire 
ropes type 3x7 and 1x7, assuming equal conditions imposed in previous 
investigations done about protection barriers. At the end is concluded that all 
models, excepting the model of internal stresses, express the physics 
present in the referred experiments and simulations thus contributing to 
augment our knowledge about the transversal impact of steel wire ropes. And 
additionally, it has been contributed with a method to determine the frontier 
conditions in a transversely impacted string that can also be used as basis to 
solve more general problems, as the oblique impact of wire ropes. 

Keywords: Transversal impact, steel wire rope, mathematical model, 
protection barrier. 
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u Displacement along the x axis 

L/j 1/2 Movement characteristic parameters of the delimiter string with 

small displacements 

Vflo Impact body velocity or its velocity at the instant of contact with 

the string or wire rope 

V Impacted string’s point velocity 

w(q) Function of a variable q 

X , y , z Coordinates of reference 

X Actual coordinate or force per unit length along the x axis in an 

ideal string 

{x,Y,z) External forces per unit of length over an ideal string 
Jo Displacement along axis y of a string’s middle or impact point 


y R.Corr.D 


Displacement while unloading of the middle point of a string with 
small displacements 



yR corr s{ed) Impact poifit displacement during the unloading expressed as a 
symmetric function of the loading of the same string with its 
unloading modulus of elasticity Ej^ and with the model of small 
displacements 

yR Con- s{ED Max) l^pac^ poInt total dlsplacsment for the string loading stage. 

yR.aprox. Impact point approximate displacement using the model of small 
displacements 

yR.corr. Corrocted displacement of the string’s impacted point with a 
delimiter movement and small displacement 

yj^j Material point displacement (model of large displacements ) 

yj^j c Material point displacement while loading (model of large 

displacements) 

y^ Material point displacement while unloading (model of large 

displacements ) 

y' \ Material point acceleration (model of large displacements ) 

y''M.c Impact body acceleration while the wire rope is loaded (model of 
large displacements) 

y' \ Impact body acceleration while the wire rope is unloaded (model 

of large displacements) 

y V Impact body velocity (model of large displacements) 

y \ c Impact body velocity while the wire rope is loaded (model of 

large displacements) 

y 'm.d Impact body velocity while the wire rope is unloaded (model of 

large displacements) 

Jb Delimiter string’s middle point displacement (model of small 

displacements) 
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Bar’s neutral axis slope value squared 

Angle between the x axis and the impacted string taking a 
triangular shape after contacted by a material point 

Finite displacement 

Infinitesimal change of area 

Infinitesimal change of the actual coordinate of a deformable 
solid or ideal string 

Small lineal deformation 

Deformation ratio 

Correction factor of the approximated displacement (model of 
little displacements) of the delimiter blue light’s middle point 

Constant 3.14159 ... 

Initial lineal density of a string 

Actual lineal density of a string 

Density [Kg/m''3] 

Stress 

Longitudinal stress beard equally by every wire in a wire rope 
Stresses parallel to the x axis along the z axis 
Stresses parallel to the y axis along the z axis 
Stresses parallel to the z axis along the z axis. 

Maximum longitudinal breaking strength of a wire rope 
Maximum breaking strength of a wire rope 



cTy M Equivalent Von Mises stress along the z axis with origin at the 

center of the contact zone between the impact body and the 
wire rope 

o'fi,Mdx Breaking strength (Pa) 

fj, Shear modulus of two steel cylinders in contact 

V Poisson ratio 



1 . Introduction. 

1.1 . Background and motivation. 

Protection barriers are used all around the world to reduce the fatal 
consequences of accidents in roads, by fall of objects in mines and 
landslides [1], [2], [3], [4], The reasons to design and install these structures 
are showed in the global annual statistics that show how more than 1.25 mill. 
Persons die in traffic collisions, that mining is responsible for 8% of all labour 
accidents and landslides caused 4626 deaths only between 2015 and 2016 
[5], [6], [7], [8], [9]. Peru, is not an exception, to the contrary, we have much 
more reasons to use protection barriers as 50% is very susceptible to 
landslides and rockfalls (with velocities up to 38.9 m/s) caused by 
earthquakes and the high frequent meteorological phenomena (El Nino), that 
caused 591 landslides and 59 deaths in the same biennium. We should also 
take into account the huge amount of road accidents: 2253 deaths in 2016 
(with velocities up to 69.4 m/s) and the local statistics in mining with 115 
accidents from falling objects in 2017 [10] [11], [12], [13], [14] [15], [16]. For 
these situations, here have been installed protection barriers made of steel 
wire ropes and steel nets in Lurigancho (Chosica, Lima), and were used 
during the summer of 2017 when El Nino caused landslides and the barriers 
stopped medium size rocks. Also, in Lima, along the cliffs of the Costa 
Verde, is installed a geo-mesh to protect from rock falls to the drivers in that 
road [17], [18]. 

In figures 1, 2, 3 y 4 are showed examples of protection barriers in 
some of its applications [19], [20], [21], [22]. 
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Figure 1: Example of a Figure 2: Road barrier 
vial barrier made of steel an impacting vehicle [20] 
wire ropes. [19]. 


Figure 3: Dynamic barrier used at ravine 
Carosio in Chosica, Lima. This stopped 
detritus during the landslides of 2017 
summer [21]. 


Figure 4: Example of a barrier stopping a falling rock 
that could impact cars if it reach a road down hill 
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However all these actions are not enough to solve our situation 
because We should not limit ourselves to wait and use these solutions, We 
must do better and new designs. But to design and improve, the behavior of 
the barrier’s structural elements as the wire ropes have to be studied as 
these elements are the principal part of many types of safety barriers as: 
flexible road barriers, rockfall barriers, and artificial bottoms used at deep 
shaft mines. In all these designs the wire ropes are impacted obliquely and 
transversely by vehicles, rocks o machinery which in turn are redirected or 
stopped to avoid or lessen economic and mortal damages [1], [2], [3], [4], 
[23],[24], [25], [26], [27], [28],[29], [30]. Due to this reality, this job, is centered 
in studying the transverse impact in steel wire ropes so, this way. We can be 
able to give more and better solutions as mentioned above. 

1.2 . Previous works. 

The development of protection barriers is based on investigations like 
the one of C. Stolle who experimented and simulated with precision a 3x7 
wire rope under transverse impact to be used in the design of flexible road 
barriers, but his numerical study is limited to this wire rope and low velocities 
(10.57 m/s) [31]. About rockfall protection barriers, Koo et al, studied 
numerically two types of this structures to predict its behavior when a 4000 
Kg sphere impacts it at 36.6 m/s, they validated its study with experimental 
data although some simulated data differed from the measured on wire ropes 
and columns [32]. Similar works were done related to artificial bottoms for 
deep shaft mines, for example, Tytko et al and P. Gospodarczyk simulated a 
design with 1x7 wire ropes in a criss-cross pattern to form a net which was 
supposed to stop machinery falling from 800m high. They were limited by the 
little knowledge about this topic and the resources to make experimentation 
at scale 1:1, hence, first they decided to test the elements or materials under 
impact individually then all this data was used to simulate and optimize their 
design [29], [33]. And the more recent work about structures with wire ropes 
under transverse impact is the one of Dqbek et al, their principal objective 
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was also to design a barrier for deep shaft mines but with a fall of 1000m and 
using 6x36 wire ropes in a criss-cross pattern testing in advance its 
properties individually in laboratory. Its main conclusion was that wire ropes 
only absorb a little portion of the energy and do not optimize its design [4], 

At this point, what can be observed in the referenced researches is how 
are limited to the designs and simulated wire ropes using only a commercial 
software LS-DYNA. Neither using any analytical model nor doing a reference 
to theoretical developments about wire ropes under impact that were being 
made by many years as G. W. Housner did in 1945 when he had to model 
the transverse impact of airplanes over wire ropes attached to balloons which 
were used as a barrier to protect a nazi installation during WWII, His model 
predicted that the airplane wings (assuming a constant velocity) could break 
the wire ropes without any mayor damage and was what happened in 
practice [34], During and after the war K. A. Rakhmatulin, also studied 
theoretically the same situation with barrage balloons, and hence, the 
transverse and oblique impact of wire ropes including friction, keeping the 
impact body’s velocity constant or with a small deceleration and assuming 
the wire rope as a string with no bending resistance and small displacements 

[35] . Then, other researches made similar formulations: F. O Ringleb, 
modeled and tested experimentally 6x19 wire ropes with a hemp core and 
used this to design aircraft arresting gear systems used in aircraft carriers 

[36] . But was N. Cristescu, who give more contributions about the impact in 
general applied to fibers and considers a broad bibliography about this field; 
furthermore he presents a numerical method to solve the case of the 
transverse impact over a string with a fixed end but assuming beforehand a 
frontier condition with a sinusoidal function which models the movement of 
the impact body (material point) [37]. Then, more recently in the nineties, J. 
Lombard, based its work on G. W. Housner’s and investigated accidents 
involving military aircraft accidents impacting transversely wire ropes of 
ropeways. He establish a purely plastic deformation and assumes a parabolic 
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profile for the string instead of the triangular one used before and also 
assumes the body velocity as constant; his results coincided with what G. W. 
Housner found and either did any contrast with experimental data [38], [39]. 
Then, what characterizes these works are the simplifications made in order to 
found an analytical solution after reducing the general differential equations 
by assuming a constant velocity, or giving a function for the kinematics of the 
for the impact body, or establishing a shape for the impacted wire rope as is 
shown in figure 5. 
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Figure 5: Idealized steel wire rope taking a triangular shape 
during the impact of a body with constant velocity v _ kv , is 
the velocity of the transverse wave generated from the impact 
point and c , is the velocity of the longitudinal wave [34], 

More than 60 years had to go by to have a formulation where is not 
assumed a shape for the string or a constant velocity for the impact body. 
That was X. Teng (2005), who modeled fracture on strings like bodies 
transversely impacted by a parallelepiped, assuming perfectly plastic 
materials and considering a high impact velocity (->500 m/s)with which 
simplifies the general dynamic equations of the string to the conservation of 
energy and conservation of momentum and founded a solution but, again, do 
not gives any experimental validation [40]. 

This way in our present We have that every experimental o numerical 
study about protection barriers containing steel wire ropes limits its findings 
to the types of wire ropes used and do not use any theoretical approach 
made before, and these ones are limited to keep the impact body’s velocity 
constant, assume an small deceleration,suppose a function for its movement. 
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or even assume a priori a shape for the impacted wire rope and constant 
properties for its constitutive material (steel). However, all this simplifications 
do not happen in real applications because the principal function of a 
protection barrier is to stop some vehicle in an accident, or a falling rock o 
machinery, and thus, preventing them from reaching people in roads, towns, 
workplaces, etc. But, this implies a great deceleration of the impact body and 
the barrier protection structure should keep its structural integrity after a 
collision(s) according to its design. 

Hence, there exists the need to design more and better protection 
barriers. And knowing the limitations that simulations, experiments and 
theoretical approaches done about impacted wire ropes, this work, is 
centered in modeling mathematically the transverse impact of wire ropes by 
any mass that could be an automobile, some machinery or a rock and this 
way it will be possible to predict, mainly: changing longitudinal stresses as a 
result of the acceleration of the impact body assuming an elastic wire rope 
which do not break or fracture in concordance with the principal function of a 
protection barrier which is protect. 

Now, accordingly to the aforementioned, a research problem is stated: 

1.3 . Problem statement. 

Develop a mathematical model to determine the emerging stresses in a 
steel wire rope transversely impacted by a punctual mass. 

1.4. Justifications. 

Justifications to work on the stated problem are: 

1.4.1 . Scientific justification. 

Because it is necessary to understand better and more the 
phenomenon of transverse impact of wire ropes. As an open field to research 
and as analytic alternative for nowadays numerical models. 
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1.4.2 . Technological justification. 

The need to predict the stresses in a steel wire rope transversely 
impacted by a mass that could be a vehicle, rock or machinery hitting some 
road barrier, rockfall barrier or artificial bottom mine, and because an 
analytical tool is needed for prediction and selection of wire ropes to be used 
in this barriers. 

1.4.3 . Social and economic Justifications. 

There is also a need for more and better alternatives to diminish or 
eliminate mortal and monetary effects due to landslides, rockfalls, road 
accidents and by falling objects in mines. Here is also included the 
developing of this solutions; so, time, financial and human resources can be 
saved when experiments and simulations will be done with prototypes of wire 
ropes under transverse impact. 

1.5 . Objectives 

1.5.1 . General objectives 

The prime goal is to develop a mathematical model for the longitudinal 
stresses, another of the the kinematics of the body impacting in its middle 
point and a last one of the internal stresses around the contact zone of a 
transversely impacted wire rope. 

1.5.2 . Specific objectives. 

After establishing and solving the mathematical models, it is sought to 
do the following: 

• Analyze the impact body’s acceleration dependence on the wire rope 
characteristics and its impact velocity. 
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• study the relation of the wire rope’s points displacements behavior 
during the impact with its properties and the impact body’s mass and 
velocity. 

• Analyze the wire rope longitudinal stresses during the loading process 
(impact body deceleration) or unloading process (impact body return 
acceleration) 

• Validate the longitudinal stresses model by comparing its results with 
the experimental data of C. Stolle’s^ research about 3x7 wire ropes 
impacted transversely and also with the numerical results of Tytko et al 
and P. Gospodarczyk’s research [29], [31], [33]. 

• Calculate the contact stresses between the wire rope and the impact 
body for each C. Stolle’s experiment and the simulation of Tytko et al, 
then, make a qualitative comparison with what was observed [31], [29]. 

1.6 . Thesis scope. 

Here, mathematical models are developed for the mechanical 
phenomena happening in a transversely impacted wire rope by a punctual 
mass. This research is limited to the impact in the point of the wire rope, 
which has an elastic behavior, with its geometry idealized as an 
unidimensional deformable string without bending strength but including a 
permanent deformation after loading to consider hysteresis. The impact body 
is represented as a material point, this is,no other effect is considered more 
than its acceleration on a straight trajectory in contact with a non-friction wire 
rope surface. This analytic approach comprises formulation of models 
oriented to predict the measures taken in experiments and simulations, i.e. a 
model offerees in the wire rope ends, a model of the impact body kinematics 
and one of the internal stresses around the contact zone. Then, validation is 
done for each model idealizing the experiments a simulations of 3 researches 

1 C. Stolle usa la nomenclatura DBC-1,2,3... para denominar a sus experimentos con cables 
impactados transversalmente y DTC-1,2,3.. para los experimentos de impaeto longitudinal [31]. 
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where 1x7 and 3x7 wire rope types were used under transverse impact 
showing a good agreement, except for the maximum internal stresses around 
the contact zone. To end this work conclusions, recommendations and 
annexes are given. 

1.7 . Outline of the thesis. 

This thesis consist in eight chapters which arrangement is as follows: 

Chapter 1 critically and succinctly reviews the social problematic 
involved and previous studies done about protection barriers. There is also 
the problem statement, justifications, objectives, the scope and a thesis 
outline. 

Chapter 2 provides a literature review including basic concepts about 
impact of solids, displacements and deformations in strings, longitudinal 
stress waves, dynamics of elastic strings, behavior of steel wire ropes under 
axial stresses and a physical description of transverse impact of elastic 
strings. 

Then chapter 3, establish the materials and methods to be used and 
followed, also gives which variables are considered along with geometric and 
physical restrictions and simplifications. Formulation and validation 
procedures are given as well. 

Chapter 4, developes the matehmatical models. 

In chapter 5, the models are validated using experimental and 
simulated data from three previous works about protection barriers which use 
3x7 and 1x7 steel wire ropes. Next, chapter VI concludes the thesis and 
gives recommendations. In Chapter VII are the references and finally chapter 
VIII includes two annexes. 
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2 . Theoretical background. 

Relevant topics for this work are: impact of deformable solids; 
deformation, displacements, stress waves and dynamics of bars and strings 
under external stresses, also,the behavior and properties of steel wire ropes. 

2.1 . <|,What is considered as impact? 

Before studying wire ropes under transverse impact is necessary some 
approaches about the impact phenomena in general, as it is not common in 
an undergraduate level. 

Impact between deformable solids is observed when discontinuities in 
properties as density or other thermodynamic quantity are present due to 
large changes in pressure by some sudden exchange of momentum and that 
then are propagated throughout the solids. These discontinues changes can 
be produced by simple percussions and we perceive them with our senses 
as abrupt and violent. Although, there is not an absolute definition of what 
“discontinuity” or “sudden event” signifies in terms of time; thus, any concept 
would be something subjective. For example, in the intellectual world of 
impact physics, a common discontinuity would be a wave front traveling from 
the contact zone to every point inside a solid that could be a bar colliding with 
another one, provoking high peaks of pressure in just some tenths of a 
nanosecond. Another example can be observed in geology, where 
earthquakes can produce severe and non-reversible changes in many 
seconds. Curiously, these cases, even happening with time and space scales 
too different, deal with the same orders of magnitude when are related with 
its life spans which are taken as its characteristics magnitudes: suppose that 
a solid medium has a characteristic space magnitud, l ^ (a bar length for 
example) and a characteristic elastic stress wave velocity of propagation, or 
sound velocity. If this body is subjected to an external dynamic load with a 
characteristic time which is, for example, the time needed to reach its 
highest level or the time when its value becomes null again; and if t^^LIc , 
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then, the stress distribution and deformation in that solid are not uniform. 
Hence, wave propagation effects should be considered. Some examples of 
characteristic times are 10 nanoseconds that is a 4.2 x 10'^® part of a human 
life of 75 years (or it could be a generation or a structure’s lifespan) while 1 
second could be a characteristic time, as it is a 6.9 x 10'^® part of earth’s life 
(or can be compared with a geological like Cretaceous or Mezosoic). These 
very little times are also why many models about solids impacting strings 
assume that after the first contact the string points change its velocity 
instantaneously from zero to the body’s velocity and is represented with a 
Heaviside function [36], [40], [41]. However, all the above only tells that every 
initial contact between solids generates a propagating stress wave and if its 
characteristic time is small hence, it is considered a percussion^. But is not 
possible to determine if those stresses are big enough to cause considerable 
structural changes inside the impacting solids. For this reason another 
criterion is considered, “The Weihrauch characterization”, there is suggested 
some references to determine how a contact between solids is clasified, 
based on its deformation ratios and relative impact velocity. If the observed 
deformation ratio has 10° order of magnitude and a velocity of 50 m/s, then 
only a localized plastic behavior is could happen, but with velocities from 50 
to 500 m/s and deformation ratios with an order of magnitude of 10^ then only 
a plastic behavior is expected. Of course these parameters are only 
referential as them do not consider many other factors like geometry, 
generated heat, friction or strength similarity between the colliding solids 
materials and its behavior with high deformation ratios as happen when a 
hard steel impacts an 1100 aluminum, this latter will develop a purely plastic 
flow typical of hyper-velocity impacts while the steel will have at most, an 
elastic disturbance [42], [43]. Hence, it is imperative to analyze every case 
and to establish which parameters are relevant. An example is how N. 
Cristescu, pointed out that if a nylon strings under transverse impact is hit by 
a body with a velocity exceeding 1/20 of the transverse wave velocity or 

2 Synonym of impact or collision, at least in mechanics. 
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deformation ratios of 33.33 m/m/s then the effect of this wave propagation or 
the local inertia forces must be taken into account. This is to say, here one 
can not simplify the string movement equations with a priori suppositions like 
assuming a triangular shape during the movement as other have done [37], 
[34], [38]. Although until now there is not a similar criterion for steel wire 
ropes that can tell when local inertia is relevant according to a impact 
velocity, this remains relatively unexplored. 

2.2 . Strings displacements and deformations. 

It is also necessary to clarify the displacements and deformation 
concepts relative to this kind of structural elements because a wire rope can 
have large displacements with small deformations. About this fact A. E. H. 
Love indicates how a bar, and by extension strings, can have relative not 
small displacements between its parts and even with this its deformations 
can be small enough to satisfy requirements of a mathematical theory in 
Strength of Materials or linear Elasticity [44]. Thus, there is an approximation 
using “small deformations” which lets consider a linear elastic response for a 
material and in the other side there is the “small displacements” 
approximation where elemental volume rotations can be considered as rigid 
body rotations and be irrelevant to calculate elastic deformations, but if taken 
into account they introduce non-linearities in our equations [45]. This latter 
condition can usually be found in conventional engineering applications, the 
classic example is the simplification of the bending moment - curvature 
relation after being linearized by considering the square of its slope [dy/dxf 
too small, compared to 1, t but this can only be done after considering small 
deflections and an almost straight beam or bar [46]. Thus, in any other case 
involving slender bars or wires where its deflections are comparable with its 
own length an elemental linear theory can not be used. Hence, to model 
bars, wires and strings shall be known that large displacements will be 
present. Which means non-linearities due to its geometry but with a material 
constitutive law that can be used to restrict the study to small deformations. 
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The foregoing have to be complemented by expressingmathematically 
the deformation of an uni-dimensional body which is an string as is the elastic 
solid studied here. There are various deformation definitions to choose from, 
but for a small deformations regime (in case of elastic metals) it is possible to 
probe that every one ends to be equal to Cauchy’s definition [47], It have to 
be clear too, that even if an impact involve dynamic stresses this does not 
necessary implies the material behavior would differ from its statics 
response. Von Karman, pointed out this on its study of longitudinal impact 
reaching plastic stresses in bars. He noticed how the ratio of deformation can 
influence stresses only at very high impact velocities, but as does not specify 
how high then the Weihrauch characterization can be used to establish a 
reference [43], [48], [49] . 


2.3 . Unidimensional stress waves. 

These appear in solids even if an interaction is not considered an 
impact. It is only needed a relative change in position between its points. 

Aiming to find a mathematical expression, figure 6, represents the 
model of prismatic bars, or straight bars, there is a bar with variable 
transverse area along its axial axis aligned with the coordinate reference axis 
X, with origin O [50], [51]. Then with an isocoric process and in dynamic 
equilibrium the following can be written: 


A(1 e) Aq , AgP^^ [EA ) „) 


( 1 ) 


d^t d^x' d^x l+sd^x 

In equation (1), a change of the transverse section area is considered 
related to the deformation in the actual bar configuration, but as indicated by 
Timoshenko, this does not take into account lateral inertia from particles 
moving due to contraction of these sections, hence,this relation is valid only if 
the stress wave amplitude is large compared with the transverse bar 
dimensions [52]. This equation can be done simpler by considering a solid 
bar with an uniform transverse section, as follows: 
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Being its general solution the one of D’Alambert: 


u{x,t) = fi{x-ct)+f2{x+ct) 


( 3 ) 
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Figure 6: Dynamic equiibrium of a deformable bar with varibale 
transverse section, x , is the coordinate of reference, X , is the 
actual coordinate; dX is a longitudinal infinitesimal element, mm ' 
and nn ' are sections, u is the displacement, 8 A is the infinitesimal 
area’s change, o is the evenly distributed stress in the cross section 
and c is the longitudinal stress wave velocity [50], [51]. 


2.4. General equations of motion for an elastic 
string. 

N. Cristescu was the first to work completely and developes in deep the 
general equations of motion for an extensible string moving in a tri¬ 
dimensional string [37], [53]. He used a curvilinear coordinate S along the 
string with an arbitrary origin, as shown in figure 7, this is an Eulerian 
coordinate and called “actual” as it measure the deformed string in the actual 
time t. The total force in any string point is denoted with T , (it is force and 
not stress because an unidimensional string has not area), is tangent along 
the string and also varies continuously with time. There is also considered a 


14 



segment between s^ and S 2 tensioned and with external forces {x,Y,z), 
which are proportional to the segment length. 


i 

/ r(sj 



Figure 7: Infinitesinmal string segment in 
a tridimensional space under internal and 
external forces [37], 

Then, the dynamic equilibrium is expressed as follows; 


*2 

f ^ 

J a,s 

Si , 

With two similar equations for y and z . Here, p is the actual string 
density and function of s and t. Projections of T over the three 
coordinated axis are ||r||(5sx/a,s), ||r||(5,y/5,s), ||r||(0,z/5ss), being 
||r|| its modulus and Sj, S 2 arbitrary, s can also be expressed using the 
Lagrangian or material coordinate , in the same way for x , y , z there is 
a expression like: x(s(so),0 , etc. Then: 
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, idem for y, z. 


Now, deformation is defined following Cauchy’s: 
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And due to mass conservation: pd^ s=pf,d^ => po=p{l+s) here, p^ 

, is the non-streched string density. Hence, the movement equations (4) can 
be expressed as follows: 


d,X , . ■ . r 

-Pa —^+(l+e)X=0 , idem for y, z. (7) 

d,t 

Here, a material’s constitutive law can be added. This is usually an 
increasing monotone function and defined as ||r||=||r||(£) . Then using (6) 
and equations for x , y and z from (7) the next system is found: 


d,s 


|r|| 

l+£ d, Sn 


a*/ d]z 
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( 8 ) 


+(l+e)X = 0 


(i+£)a,^s„ ^"ar''a,^s„ (i+e)' a^^s„ 

With similar equations for y and z . This partial differential system of 
equations is hyperbolic, of order 2 and has 4 different families of 
characteristic curves, which are: 
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d,t — 

Po(l + f)'^ 

’ d,t - 

Po d^e 


1/2 


( 9 ) 


This is why characteristic lines are generally curves with a slope as a 
function of ||r|| and e . From a mechanical point of view this lines represent 
wave fronts, then there are two propagating waves with the following front 
velocities. 




Po(l+e)' 
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1 d, 


Po d„£ 


1/2 


( 10 ) 


Here is the transverse wave’s front velocity and only change the 
string shape but does not deforms it, that is, a material point reached by a 
transverse wave acquire only a velocity to its tangent. But if a point is 
reached by a longitudinal wave front with velocity c„ it cause only stretching 
but not modify its slope. 
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2.5 . Modulus of elasticity of steei wire ropes 

As the foregoing theory is used to model steel wire ropes, it is 
necessary to know how this mechanical elements behave. Elongation in a 
wire rope depends, of course, on how its parts are arranged and what 
materials are used (Figure 8) specially its wires, because, if its core is 
synthetic then it is not taken in count as does not contribute greatly to its 


elongation properties. 



Figure 8: Steel wire rope 
elements [54]. 


Generally speaking, the elasticity modulus of a wire rope is very 
different from its materials modulus. They stress^-strain curve is usually non¬ 
linear, this is, given a type of wire rope its modulus would not be constant 
above certain stress and is a function of the level of tension. This is why in 
practice is used an average elasticity modulus. Although, simple spiral ropes 
made only by twisted wires have a minimum non-linearity and can be 
neglected in most cases [54] [55]. For simple spirals and some wire ropes 


3 Stresses in a wire rope are defined as o=T/A, here A is the metallic area or a combination of 
transverse areas of every wire. 
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made of many strands with rounded wires it is possible to approximate its 
modulus of elasticity using analytic methods only, but this is not too exact for 
most wire ropes with more elaborated designs and experimental data is 
needed in order to approximate its modulus that depends on what level of 
load and how was loaded before. Additionally, the elastic limit can be defined 
as that load from where a stabilized wire rope cannot reach its previous 
length after unloading. Then, some definitions of the elastic modulus which 
have practical relevance are [54] [55]: 

• £s(cr,o„er> <7upper) - fhls Is the slope of a secant line between an initial 
and final stress. It is used with fluctuating loads. An special case 
occurs when the initial stress is zero Es(o,a'„pper) ■ 

• Et^^,p{a),Et^aown{o‘), these are the slopes of tangents in any point of a 
loading and unloading curves respectively. 

Its important to know that a calculated elasticity modulus does not have 
practical meaning if previously the wire rope has not been unloaded to the 
load from where it will be used. This has to be done because wires rearrange 
inside and also friction takes place with every loading and unloading, not only 
elastic deformation occurs. For these reasons it is recommended to stabilize 
a new wire rope, before making any measure, with 10 loading cycles. It is 
also possible to use the tangent slope in every point but in practice is only 
used as an evaluation parameter to make comparisons. Figure 9 shows the 
geometric interpretation of every modulus definition, slope of line OB is 
modulus Es(o,<7„pper) between a zero stress and the maximum at which the 
wire rope works. If there is an initial stress (point E) and a final one (point F), 
then the elastic modulus would be the slope of the secant line EF and equal 
to E^laE^o’p) ■ ■ Also, in the same figure, tangent lines Ot and in point A 
Ef o{a),E^ . And here a correct or practical modulus would be that 
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measured as the slope of line CB where the wire rope has been unloaded 
[54], [55], 



Figure 9: Elasticity modulus definitions as the slope of 
secant lines of OB, CB, and EF, and of tangents in O or A. 
a, is the stress, e^, is the final deformation, s^, is the 
deformation while loading y , is the deformation while 
unloading cp , are the angles of the tangent and secant 
lines. [55]. 

Next, Figure 10 illustrates, in a general way, the non-linear behavior of 
many new steel wire ropes when being subjected to one cycle of loading and 
unloading. Plastic coated wire ropes or with a fiber core tend to have a 
completely non-linear behavior while spirals (strands) have a predominantly 
linear response. An important observation here is how after a first cycle every 
type of wire rope presents an increment of its elasticity modulus, this is why 
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10 cycles are recommended before an stable measure is taken. In the same 
figure the lineal idealization of a loading and unloading cycle is represented 


fora new wire rope using straight lines [31], [54], [56], [57], [58], [59], [60]. 



steel bar not drawed. 

Linear idealization of a 
steel wire rope. 

Wire ropes of types 1x, 
3x. (with only 1 strand 
or 3). 

Wire ropes with more 
than 3 strands: 4x, 6x, 
8x,..18x..., filled with 
plastic or with a fiber 
core. 


Figure 10: Curves stress o - deformation s , of some types of new steel 
wire ropes, after being loaded and loading for the first time. [31], [54], [56], 
[57], [58], [59], [60]. 


20 







2.6. Physics of elastic strings under transverse 
impact. 

Here, only geometry and English are used, avoiding mathematics on 
purpose, to communicate the involved physics happening when a body 
transversely impacts a wire rope, and also how its quasi-isomorphic 
idealization, that is, an uni-dimensional string without bending strength 
behaves after being impacted by a material point. 

First, physics of a transversely impacted steel wire rope can be 
describe as follows: 

1. When a body impacts transversely a wire rope then the initial contact 
is produced over some or many superficial wires which in turn will 
deform and transmit the contact force all along the wire rope. From 
this first instant there is also a relative velocity between the body’s and 
wire rope’s contacting points. This contact also means there surge 
stresses, shear and compression, apart from the longitudinal ones. 

2. The longitudinal deformation appearing in each wire correspond to the 
one of a double or simple helix body (as in a spring), depending of 
what type or location it has. This means, there is a tension, torsion, 
bending and shear stresses [54]. 

3. Because of the foregoing, surging stress waves are also of torsion, 
bending and also of shear nature. Every one of these perturbations 
travels through the curve axis of every wire and they add up each 
other to produce the wire rope’s total displacement. In addition, over 
the surface of every wire appear shear waves due to friction. And all 
them travels up to each wire rope end where refraction and reflection 
occurs. Then every reflected wave goes back and combines with the 
actual stresses in every wire changing continuously as the contact 
force or reaction force over the impact body changes. 


21 



4. Intern friction forces between wires are not a function of a constant 
coefficient only but an state function of each material in contact [61]. 

5. Also, medium resistance and gravity should be taken in count as them 
influence the wire rope geometry. 

Hence, if any research considers not to simplify the wire rope geometry, 
the model should include all deformations with its internal and external forces 
mentioned before and how reflection and refraction happens in each end. 
Interaction with the impact body should be considered too, as its dimensions 
can be comparable to the wire rope length and because friction arises in the 
contact zone. But this way is not at the author’s reach, by now, therefore here 
a wire rope is modeled as an uni-dimensional string and the impact body as a 
material point. 

Thus, in order to idealize the wire rope movement these general 
simplifications are made: 

1. The wire rope is reduced to its geometrical axis which initial state is an 
uni-dimensional curve, straight with uniform density resting in a plane. 

2. Only the axial deformation is considered, choosing the Hooke’s law to 
relate it with the axial stresses. 

3. The impact body is reduced to a material point with a given velocity 
which starts to decelerate when it touches the wire rope. 

4. It is considered that the point where the body impacts, acquires the 
same velocity instantaneously after the first contact. 

This last simplifications is a consequence of assuming a deformable 
and with no bending strength string. It is a discontinuous change of velocity 
which also means the geometric model has a singular point, kink or 
discontinuity on its tangent at the impact point. This singularity is where a 
transverse wave front begins provoking that all contiguous points gradually 
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get a transverse movement which means a continuous change of the total 
mass in movement. Its velocity is represented with the symbol Cj of 
equation (10). 

In the other hand, propagation of stretch is done longitudinally and its 
velocity depends on how much deformable is the string. Here depends on its 
elasticity only and represented by Cjj of equation (10); it separates relatively 
every point it reachs but without altering its tangent direction as the 
transverse wave does. 


Figure 11 shows how the shape of an idealized wire rope changes 
(before any reflection takes place) after being impacted transversely by a 
punctual mass. The transverse wave front has a velocity Cj which varies 
with tension, The longitudinal wave front velocity depends on how tension 
changes with respect to deformation, this is, in the elastic case it is a 
constant Cj^. The impact body’s initial velocity, at the instant of impact, is 
Voo, the body has a mass M , then V is its actual velocity, , is the initial 
coordinate (Lagrangian) for every string point and S is the actual one. The 
space where the movement occur does not offer any resistance or friction. 


M * ^00 Transversal wave front 




Longitudinal wave front 



Wire rope idealization as an 
^ V Unidimensional string 


Figure 11: Movement of a string transversely impacted by a punctual 
mass. M, is the impact body mass, Voo, is its initial velocity, V, is its actual 
velocity, O, represents the origin of the coordinated system x, y, S, is the 
actual curvilinear coordinate and So, is the referential coordinate. 
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Ideally, physics of the transverse impact of a material point with a 
deformable unidimensional string has the following stages: 

1. Stage 1: At time zero an initial contact is produced between the 
material point and an point in the string, and instantaneously both get 
the same velocity. This way the loading stage begins. 

2. Stage 2: Here the body starts to decelerate and from the contact point 
a longitudinal wave travels with the velocity of sound along the string, 
at the same time a transverse wave, that only changes the tangent of 
every point that reachs, propagates to both ends. This last wave is 
how local inertia appears although the longitudinal wave also 
contributes but in minor scale because it stretches the string and 
changes its local density. 

3. Stage 3: If the body continuously decelerates and the longitudinal 
wave fronts reach the ends then the reflected waves are added to the 
stress state that also changes with the reaction force over the impact 
body, this is, the stress wave produced initially in the contact point 
comes back and influence the movement. The transverse wave also 
reflects or refracts and begins to influence the string points and body 
movement as it interacts with the wire ropes inertia. 

4. Stage 4: The wire rope gets a zero velocity and begins to change the 
sense of its movement (keeps its straight direction). This way the 
loading stage finish and the begins to unload. Then, should be 
considered if the impact body is kept attached to the wire rope and if 
this one has a permanent deformation, at some time while returning, 
the point of contact will have a minor velocity than the material point 
and its movement will continue until its energy totally dissipates. Here 
the movement is considered only until the wire rope tension becomes 
zero again indicating the end of the unloading stage. 
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3 . Methods and materials. 

3.1 . Information gathering. 

All the necessary information to establish the problem and solve the 
equations is found, mainly, in works done by F. O. Ringleb, N. Cristescu and 
K. Feyrer [36], [37], [53], [54], Validation with experimental and simulated 
data is done by using C. Stolle, Tytko et al and P. Gospodarczyk’s 
researches [31], [29], [33]. 

3.2 . Object of study. 

Here is studied a wire rope anchored in both ends, made totally of steel, 
helically braided, and with circular cross section wires. It gets a movement 
after being impacted transversely by a punctual mass. And in order to 
analyze mathematically this phenomenon its geometry and movement are 
idealized applying many restrictions and simplifications. To validate the 
established mathematical models data from C. Stolle, Tytko et al and P. 
Gospodarczyk researches about transversely impacted wire ropes type 3x7 
(ASTM A741 -11 (2016)) and type 1x7 (BS 5896-2012) are used [31], [29], 
[33]. 

3.3. Methodology to establish the mathematical 
models. 

The method used in this research is part of the mathematical modeling 
methodologies and these in turn part of general modeling which main 
purpose is to relate the real world with some simplified physical and 
conceptual representation moving back and forth cyclically between reality 
and imagination. Specifically this work is within the Philosophy of change 
founded by Newton and developed, here, by J. Olivencia'*. 

The first step to establish a mathematical model is to set restrictions 
and simplifications based on the real phenomenon and then mathematics 
can be used as language to express what is in the idealized world [62]. In this 

4 Professor of mathematics at the Universidad Nacional de Trujillo, Peru. 
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first part, the idealization is done following to P. Teodorescu. Thus, 
observations of transverse impact in wire ropes in many types of protection 
barriers are translated to an ideal reality in the mechanics of deformable 
solids modeling. This way, any physical bodies using its more simplified 
shape: bars are reduced to its axis, plates are reduced to its middle surfaces, 
or a wire rope is taken as its longitudinal axis. Hence, all is framed into the 
Strength of Materials [63]. Then at the ideal world and needing to 
establishing a theory, the Rational Mechanics paradigm can be followed to 
set specific steps to model mathematically a wire rope impacted transversely 
by a punctual mass. So, following what A. Klarbring pointed out the method 
is: [62]: 

1. First, the phenomenon is idealized setting geometrical and 
physical restrictions and simplifications. In this way the bodies 
under study are described with a geometrical and mathematical 
language. There are four models for material bodies: discrete, 
unidimensional, bi-dimensional and three-dimensional ones. 

2. Second, the universal laws should be included: mass, energy, 
linear and angular momentum conservation; then, expressions in 
terms of differential equations can be written accordingly to the 
bodies under study. 

3. And the third step to obtain a solvable mathematical expression is 
to introduce particular laws: physical and mathematical 
restrictions, force and constitutive laws (that can be experimental 
or idealized) 

Furthermore, a validation methodology implies a quantitative and 
qualitative comparison between the experimental or simulated results and 
what is predicted with the developed theory. And then analyze why there are 
differences and coincidences, taking in count how the chosen parameters for 
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this model, the wire rope and impact body characteristics influence the 
analytical and experimental results. 

3.4 . Research variables. 

If the mathematical models of this project are classified as described by 
K. Velten y a F. E. Cellier, the system in study in each model is: technical, 
deterministic, continuous, bi-dimensional and physical. Knowing what has to 
be answered each model is mechanistic, instationary (in time), distributed (in 
space), inverse, of research (not management), analytical and in macro¬ 
mechanics. And the mathematics is: not-lineal, numerical, continuous, and 
with differential equations [64], [65]. In sum, all models are gray-box models 
as represented in figure 12. Input variables are: 

1. Mass of the impact material point, . 

2. Material point velocity at the instant of impact, . 

3. Wire rope’s constitutive law: Hooke’s law for the loading and 
unloading stage, a=Es. 

4. Wire rope Poisson ratio, v. 

5. Wire rope transverse or metallic area, A . 

6. Wire rope equivalent diameter, D . 

7. Wire rope lineal density, p. 

8. Wire rope straight length: L . 

The controllable variables are: 

1. Wire rope elastic limit. 

2. Wire rope breaking strength, ■ 
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Uncontrollable variables are:. 

1. Margin of error of the solution method. 

Output variables are: 

1. Wire rope stress or force, a,T . 

2. Acceleration of the material point, y' V ■ 

3. Velocity of the material point,, y . 

4. Displacement of the material point, . 

5. Equivalent stresses around the impact zone the wire rope, . 


CONTROLABLE VAEUABLES 




MATHEMATICAL MODEL 
GEOMETRY 

GENERAL LAWS OE MECHANICS 

CONSTITUTIVE LAWS 

NUMERICAL AND ANALYTICAL SOLUTION 




UNCONTROLLABLE VARIABLES 


Figure 12: Grey box representation of the established 
mathematical models. 
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OUTPUT VARIABLES 


3.5 . Procedure. 

The procedure follows the outlined method before. This way, first 
geometric and physical restrictions and simplifications are established. Then 
a procedure is formulated to establish the mathematical models, its solution 
and the validation with experimental and simulated data. Physical and 
geometric simplifications and restrictions. 

This study has as reference the first impact scale from the Weihrauch 
characterization, this means plastic deformations are localized (contact 
zone), impact velocities are under 50 m/s and deformation ratios with an 
order of magnitude of 10° mm/mm/seg are expected [43]. Restrictions and 
simplifications for the steel wire rope are the following: 

1. The wire rope geometry is simplified to an unidimensional string with 
zero initial tension and uniform density. A transverse area is assigned 
as a constant property in every point. 

2. Only wire ropes made of steel are modeled, with helically braided 
wires and a circular section and without natural or plastic core. 

3. Only linear effects from the stress-strain curve are taken in count and 
up to the wire rope breaking strength 

4. Effects of torsion and shear stresses in the wires are not considered. 
While its friction and wires rearrangement, during the loading and 
unloading stages, are modeled using two different elasticity modulus 
for each stage. 

5. Small deformations are supposed and only dependent of time but not 
of position, temporal ratio or temperature. Its material is considered 
isotropic. 

6. It is also supposed that the applied tension does not change the wire 
rope volume, this is, the steel suffers an isochoric process and does 
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not present any shock wave that could provoke a phase change of the 
materials considered. 

7. The impact is produced in the wire rope middle point or coordinate 
zero in the string. 

For the impact body the following is used: 

1. All bodies that could impact as: rocks, vehicles or machinery are 
assumed as material points and with a known mass and velocity, 
although in the contact stresses model it is assigned a shape to allows 
calculations. 

2. Its movement direction is always straight and transverse to the wire 
rope at the beginning. 

Restrictions for the wire rope and body movement are the following: 

1. Gravity and medium resistance (air) are not considered. 

2. The point of the wire rope were the body makes contact changes its 
velocity instantaneously from zero to the body’s velocity. 

An idealized wire rope impacted by a material point is presented in 
figure 13, mass of the material point is M , its initial or impact velocity is 
Voo , the reference coordinates are , actual coordinates are S and the 
body’s actual velocity is V .Procedure to formulate the mathematical models. 



Figure 13: Geometry of an idealized wire rope transversely impacted by a 
punctual mass. M, is the impact body mass, Voo, is its initial velocity, V, is its 
actual velocity, O, is the origin of coordinates x, y, S and So, are the string’s 
actual and referential curvilinear coordinates L, is the wire rope length. 
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1. First, a continuous time t is defined, symmetric in every moment 
(Thus, conservation of energy and mass can happen) and independent 
of the space or existing bodies in it. An space is also defined from an 
inertial reference with cartesian coordinates Oxy , euclidean 
properties, bi-dimensional, continuous and symmetric (Thus, 
conservation of momentum can happen). 

2. In the defined space is a steel wire rope idealized as an uni¬ 
dimensional finite solid, deformable, with uniform density, without 
bending resistance, of elastic material, initially straight and coincident 
with X axis. 

3. There is also a material point defined with mass M , and is the 
idealization of the solid body, and reaches the wire rope with an initial 
velocity at its middle point in a negative direction along the y axis 
and it is presupposed that will have a large deceleration at least during 
the loading stage. 

4. After defining the geometry, the general laws of mechanics are used to 
establish mathematically the general model, these are given in 2 the 
theoretical background, hence, here only the constitutive law is 
needed. This law is defined with an equation relating the stress, strain 
and the elastic modulus a=Ee or its equivalent with the string 
tension using Hooke’s law T=AEe with a modulus of elasticity if it is 
loaded and a different when it is unloaded. 

5. Next, the equations of motion are stated for the impact body, doing 
mathematical simplifications in concordance with the physical 
conditions. 

6. After finding numerical answer using the experiments and simulations 
from the references mentioned, a model is established to predict the 
stresses inside the wire rope around the contact zone, or at least the 
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maximum ones. To do this Hertz theory is used and a geometry and 
material is assigned to the impact body. 

3.5.1 . Procedure to validate the mathematical 
models. 

Here only a quantitative comparison can be made for the longitudinal 
stresses at the anchorages and the body kinematics because the 
experiments and simulations referenced only present these data. The contact 
stress model can only be compared qualitatively with what was observed 
because there are not numbers recorded for each experiment or simulation 
about the internal stresses in the wire rope around the contact zone. The 
following is considered: 

1. If the wire rope experimented or simulated is new, then its stress-strain 
curve have an initial non-linearity (friction and wires rearrangement), 
then there will be a linear part followed by a small non-linear transition 
to the plastic zone up to it breaks. Therefore here, as was established 
before, is considered only the linear part but keeping the same 
breaking strength which was measured (in the referenced works) or is 
standard. 

2. Calculations are done from the first contact or after the simulation or 
sensors detect a tension bigger than zero and up to it is zero again or 
the wire rope breaks. 

This is how, the longitudinal stresses model and the contact stress 
model are validated with the experimental results of a 3x7 wire rope under 
transverse impact done by C. Stolle and with the simulations of a 1x7 wire 
rope under transverse impact done by Tytko et al and P. Gospodarczyk [31], 
[29], [33]. 
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4 . Mathematical modeling. 

4.1 . Modeling the longitudinal stresses. 

To solve the differential equations (6) and (7) which represents the 
general movement of an uni-dimensional string, many methods have been 
developed along the past century®. The last one is that of Yu. A. Dem’yanov 
(1993) who presented an asymptotic method (without giving examples or an 
experimental validation) based on the method of characteristics to solve the 
problem of wave propagation over a transversely impacted string by a 
material point [66]. In order to find an answer he makes some simplifications 
like assuming the string deformation and the material point deformation as 
functions of time as follows. In figure 13 are the variables, Bi are constants: 

e{Sot)=Bt'" ; B,m>0 (11) 


y(So=0,r) = y„(0 = -Bif";Bi>0 (12) 

So, if there is a perfectly elastic string with small deformations (Which 
does not limit the possibility for large displacements) and equation (11) is 
used with Hooke’s law T=EAe , then, the movement equations are: 


d,s 


l|r|| 5s/ 


l+e a. So 




T d Y 
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l+£ d, s„ 


-Po 


d/ 


= 0 


(13) 



And both equations can be solved with a change of variable [67]: 

>f.T=w(q);q=|a(2+fii)'(S„+B2)'-4r''''”|''';a=^ (15) 

A common solution is given for both: 

5 Before, many others did experimental and analytical research about strings (wires) under transverse 
impact using other methods to solve the partial differential equations (6) and (7), but assuming a 
constant velocity for the impact body or as an exponential function of time. See: “Experimental 
evaluation of material Behavior in a wire under transverse impact”, A. B. Schultz, P. A. Tuschak, A. 
A. Vicario, ASME 1967, “Determination of the stress-strain relations in a wire under transverse 
impacf’ A. B. Schultz, T. C. T. Ting , ASME 1968, “Determination of the unloading boundary in 
transverse impact of an elastic-plastic string”, P. A. Tuschak, A. B. Schultz ASME 1972. 


33 




X. Thus, here is a contradiction of what its assumed by the proposed 
method and what results by using its assumptions. And even more it does not 
give a solution about how the unloading stage has to happen or how it could 
be modeled. This incoherence has its roots in how the geometry is idealized: 
from a solid to a material point and a wire rope to an uni-dimensional string 
without bending resistance. An explanation is presented in Annex 1®. 


In the other hand, the experiments and simulations, referenced here to 
validate the models, only give forces at the ends of the wire rope with the 
acceleration, velocity and displacement of the impact body. Hence, these 
frontier conditions have to be determined first, and with that the stress state 
is determined in every point and in its critical points: ends and impact point. 
In sum, the objective to know in deep the movement of a transversely 
impacted wire rope have to be delimited to what was measured and 
calculated experimentally, and in general to how the engineering application 
is made. Then, with these formulations proceed to make a validation. This 
means the specific searched values are the forces at the ends and the 
displacement, velocity and acceleration (or force over it) of the impact body. 
After validating these results the stress values inside the wire rope can be 
approximated considering what predicts the internal stresses model around 
the contact zone. But, first, to reach these answers it is necessary to observe 
specific and common physical characteristics of each experimental and 
simulated case [29], [31], [33]: 


6 N. Cristescu avoided this problem by assuming a priori a sinusoidal function to model the impact 
body displacement as a frontier condition and then found the string movement and stresses [37]. 
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1. The mass of the impact body is very large compared with the wire 
rope mass^ 


Research of: 

Wire rope 
type 

Wire rope mass 
(Kg) : M „ 

V o/ wireRope 

Impact 
body mass 
(Kg): 

Ratio 

M /M „ 

c WireRope 

C. Stolle 

3x7 

1.06E+01 

7.76E+02 

7.34E+01 

Tytko et al. 

1x7 

1.17E+00 

4.00E+03 

3.42E+03 

P. Gospodarczyk 

1x7 

8.78E+00 

1.80E+04 

2.05E+03 


Table 1: Comparison between the wire rope mass and the 
impact body mass. 

2. Time to reflection, or time taken by the longitudinal stress wave front 
to go from the impact point to an end (anchorage) rjj=0.5L/C„ , is too 
small compared with the duration of the movement lifespan or loading 
and unloading stages together. This relation is shown with the ratio 


ttotai/tR in table 2: 


Reasearch of: 

Wire 

rope 

type 

C||, Wire rope 
sound 
velocity 
(m/s). 

Wire rope 
half length 
0.5L (m). 

Calculated 
time to 
reflection t^ 
(Seg) 

Movement 
lifespan t,„,^, 
(Seg). 

Ratio U,/ 

tR 

C. Stolle (exp. 
1,2-4) 

3x7 

3.17E+03 

4.31 E+00 

1.36E-03 

3.00E-01 

2.21 E+02 

Tytko et al. 

1x7 

5.00E+03 

5.00E-01 

1 .OOE-04 

3.80E-02 

3.80E+02 

P. Gospodarczyk 

1x7 

5.00E+03 

3.75E+00 

7.50E-04 

1.00E-01 

1.33E+02 


Table 2: Time to reflection taken by the longitudinal stress waves in a transversely 
impacted wire rope compared with the wire rope movement lifespan. 

Then, in accordance with the ratios of tables 1 y 2, it is possible to say 
that in the described movements the following happens: 

1. The wire rope inertia is too small compared with the impact body 
inertia. This means that reflection of the transverse waves does not 
has a major effect neither in the forces inside the wire rope nor the 
impact body movement. 

2. As the time to reflection is too small compared with the total time of 
loading and unloading, ttotai, of the wire rope. And apparently at first 


7 For P. Gospodarczyk’s simulation the mass of every wire rope in its simulated base is 
considered, this is 420 1x7 wire ropes of 7.5m length and 0l5.7mm each one. 
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glance it could be considered that every time a time, tR, passes, the 
stress waves make the wire rope tension uniform along it from each 
end to the impact point. 

Hence, point 1 permits to establish the equations considering a 
concentrated mass of the wire rope at the impact point instead of using a 
distributed one, because its inertia does not contribute significantly to its 
movement. Point 2 says that the time taken to have almost uniform tensions 
along the string is too small that the tension in each end is very similar to the 
tension in the impact point. However, because the wire rope shape, during 
the movement, is not a perfect triangle, then, the force at the impact point 
has a different direction depending of what angle has the wire rope there. 
Hence, to model the forces at the ends and the kinematics of the impact body 
a geometrical model is established as^follows: 



Figure 14: Model of delimiter ideal strings delimiting the idealization of a 
transversely impacted wire rope on its central point. M, body mass, Voo, its 
initial velocity, V, it’s actual velocity, O, origin of axis x, y, So, reference 
lagrangian coordinate are the anchorages. Yr, displacement of the light 

blue line. Ym, displacement of the green string’s middle point (or the material 
point) and L is the distance between anchorages. 

In figure 14 is used an analog concept to that of frontier conditions but 
here are presented geometrically as “delimiter movements”. This way the 
string movement is delimited between the movement of two strings which do 
not have a distributed mass, are impacted in its middle point by the same 
body and have these characteristics: the light blue string is tangent to the 
ideal string at its ends hence tension in these points is the same for 
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both strings; the green string, which is not tangent to the ideal string in any 
point, has the same extremes and its middle point S„ = 0 coincides during the 
movement and does not have also a distributed mass. Both delimiter strings 
shapes are always triangular during the loading and unloading stages, and to 
consider the contribution of the wire rope inertia the impact mass is redefined 
as a combination of the impact body mass, , and the wire rope mass, 
^wireRope ' Fuilhermore, it is assumed that the light blue string has small 
displacements while the green one has large displacements to be able to 
model the impact body movement. This does not mean the green string 
tension approximates the tension at the impact point of the ideal string 
because its triangular shape implies that its deformation is larger at that 
point. This green string is a reference to model the movement of the impact 
body. Thus, the geometrical model of the light blue string can be written as a 
mathematical model taking as a reference the movement of a mass M 
attached to two linear springs as in figure 15. 


Left spring tension level. 



Figure 15 Analogy between the delimiter model of tension at the ends of the string and 
a model of ideal springs. T, is the springs tension, a, is the springs slope angle, F, is the 
reaction force over the impact body, M, is the mass of the impact body, Yr, is its 
displacement, x, y are the coordinates and L, is the anchorage distance. 

Idealization of linear springs does not consider its real geometry and the 
model says the force felt by the ends and mass M, T, is 

proportional to the displacement y,^ and a function of time. This is, 
transmission of force from the impact point and the string ends is 
instantaneous and furthermore the stress-strain relationship is always linear 
for the loading and unloading stage as shown in figure 16. 
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model 

with 


Figure 16: Stress (a)-strain (e) curve of a loading^his movement 

and unloading cycle of a perfectly elastic wiresmall 

rope. 

displacements is: 



F=2Tsena ; sena= 


ylL^/4+y] 


F=2 


yR _ jfi 
con ——» , 

^|L^/4+y] 


The last equation, (19), is valid only if tension T is uniform all along the 
wire rope or string. Now, with this point of reference in the impact over an 
ideal string the tension levels are as follows: 

1. Force is a function of position and time: r(So,r), with 
S„ e [-L/2,L/2]. 


2. Because anchorages are considered perfectly rigid then on this points 
the tension is duplicated (reflection in a rigid frontier). 

And here is presented a geometrical model to show how forces change in 
a delimiter string which has, in every moment, a triangular shape and a 
concentrated mass on its middle point. Changes are presented for some 
multiples of the time to reflection rjj=0.5L/C„ with Cjj=V£//o : 
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Idealized wire rope Rigid anchorages 


T(-L/2,0.5g=0 C„ T(L/2,0.5g=0 

iF(ig 



T(-L/2,1.5tJ 


T(o, i.sg T(o, i.sg 

A F(2tJ 


S T(L/2,1.5tJ 



■T(-L/2,2g 


T(0> 2g T(0, 2g 


T(-L/2,2g) 



T(-L/2, 3tJ 


T(0, 3t ) . T(0, 3tJ 


F(4tJ 


T(L/2,3g > 



► T(-L/2,4g T(0,4g C„^C„ T(0, 4g T(L/2,4g) 

tension level Overlap and add each other along the "wire rope 

Right side tension levelj 


Figure 17: Change in tension along the points of a wire rope idealized as 
a string impacted in its middle point and taking a triangular shape. M, is 
the mass of the impact body, a, angle of each right side of the string, F, 
reaction force over the impact body, T, distributed tension over the string. 
Cm, is the longitudinal stress wave velocity front velocity. The wire rope is 
idealized as a string without mass and with a “delimiter movement”. 
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According to figure 17, it is not until time that the string has 
stresses in all its points and, in this case, are symmetric with respect to the 
middle point. But it is not until time 2tg that tension can be considered 
uniform in both halves, this is, leveling of tension happens after a cycle of 
reflection when the wave front, with a zero stress level, is back at the impact 
point. Then, the same happens at time next at time 6tjj and so on. 
Thus, it can be say that from time t=2tjf and after time intervals At=2tjf, 
tension in the wire rope can be considered uniform along it. Hence, equation 
(19) is valid after the same time intervals 
2it,i ; i = i e IN-(0} . But as time to reflection is 

equivalent to travel L/2 and because the impact point coincides with the 
middle point then the uniformity happens after a distance 2Cjj{itj() . After 
this, the next expression can be written: 


F=2AE 


yn 


= 2AE 


yn 


= AE 


Tr 


2(i[L/2]) ■■■“2(iC„(,) ““c„(iO 
Another way to express this model is “the transverse movement of a 
material point attached to two springs as in figure 15 is valid only after time 
intervals At=2tjf ”, this is, each time this interval passes the next happens: 


T{-Ll2,it^)^T{0,it^) ; Vi G IN-{0} (21) 

As the time interval At=2tg is very small compared with the time for 

loading and unloading, then, it can be considered that equation (21) is valid 
for these two stages. Hence as also comprise all this time, varying i is 
possible to approximate the movement in every instant by writing and 

the following: 


F=AE^^ = AE — 


( 22 ) 


The physics expressed with geometry before can be translated to a 


natural language using terms from the method of changes used by Newton: 
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“Each change of the reaction (contact) force F between the wire rope 
and the impact body, this is, F+d,F{At) causes a change in the wire rope 
tension, which is r+d,r(Ar). Then this value travels from the impact point 
through the string with the velocity of sound C„ (for that wire rope 
construction and material), after a time it reachs and reflects at the ends, 
considered perfect rigidly here. There the wave duplicates its amplitud, 
reverses its sense of motion but not its nature (it is still tension®). This means 
the waves from each side of the string, here equals because the impact is 
over the string middle point, after a time, overlap and add when they 
coincide along the string points. But, as the time taken by each change of T 
to cover all the string is also very small, then is considered that the level of 
tension after time intervals of Ar=2tg is uniform and by extension to every 
instant t. Also, because neither internal friction nor some damping is 
considered this results in a resonant wave motion” 

In equation (22), division by Cjjt means that point 5^=0 , or the impact 
point “see” changes traveling a distance C„r and not simply L/2 . While 
division by 2 of L here is turned into multiplication by 2 of the force which, 
geometrically, is the duplicative effect of reflection at the perfectly rigid ends, 
and it is also equivalent to travel half of the string back and forth by the wave 
front or the entire length in once as the impact is at the middle point. 

Now, it is possible to apply Newton’s laws to model the movement of the 
material point. And the next differential equation is obtained: 


dfvR Yb 

M-!^ = -AE-’' 


dX 


Cijt 


(23) 


But as was observed with the use of Yu. A. Dem’yanov’s method and 
the work of N. D. Cristescu, its equations can not be integrated from time 


8 Reflection of a wave is equivalent to the overlapping and addition of two identieal waves, being of 
tension or compression, traveling to each other and when they met, the points reached by both are 
stretched or eompressed by both waves and as they are equal then the amplitude is duplicated. 
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zero. The same happens here, thus,the movement is considered from time 
[66], [37], And equation (23) becomes: 


d,t' 


AE yn 
M L + Cjjt 


(24) 


This equation, (24), is of the Emden-Fowler type and according to K. S. 
Govinder and P. G. L. Leach, this type can be integrated and has an 
standardwhich is [68] 




In equation (24), function f and coefficient n are: 


(25) 


m=- 


AE 


n=l 


M(L+C„r) ’ " " 

Then, knowing that a solution exists, an approximate solution analog to 
the lineal solution of equation (19) can be established, and should has the 
following shape: 


T B. aprox .=-63 sen (VB4 0+^5 cos (V Bfi 0 
Its initial conditions are: 


(27) 


; ^('=o)=v., (28) 

The mass to consider is . And an approximate solution 

for the displacement (model of small displacements) of the light blue string’s 
middle point is: 


Tb. 


aprox . 




EA t 
M ^L+Cjjt 


(29) 


Equation (29) describes a resonant system which amplitude grows and 
its frequency decrease continuously with time, this is, only represents a 
solution for a first loading stage. And, to test if this is an approximate solution 
it is replaced in equation (24) as follows: 
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d^tyn 


.aprox. _ 


u; yn. 


aprox. 




aprox 


dX 


4 L + Cjjt 


L+C,jt 




ui 


4(L + C„0 


9{t) 


; U,= 2 M 

^ 2 y EA ' y M 


( 30 ) 


Where: 


i(0= 


L + Cjjt 


4Ut 


C„t 



+ 

2 

4 

L + Cjjt 1 


(31) 


4AE U,U,C],t _ 

— '/ - z9(t J = — 

M{L+C,jt) 4{L+CjjtY 


U2t 


Vl+Cjj< 


(32) 


yR^yn . aprox . j{t)+9{t) 


(33) 


9{t)=- 


[/^(L+C„0 


cos 


U^t 


(34) 


^JL + C,It 

If every constants are evaluated in j{t) and g{t) for r<l then: 
L^M/Ap, as E has an order of 10“ [ Pa], M and order of 10^ [Kg], A 
of 10“''[fii^] , p of 10^[Kg/rn], C„ of 10^ and V up to and order of 10^. 
And with all these: 


j(0>i ; 


4AE 




(35) 


M{L + Cnt) 

And j{t) , with the experimental dimensions and suppositions its value 


tend to be 1.1091. Hence a correction factor can be established for yR^^prox 


9 . 


L 

J j{t)dtt 


(36) 


'Corr. 


= lim 


= 1.1091 


9 This calculation can be probed using QtiPlot and the attached files. 
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Thus, the corrected displacement of the string’s impacted point with a 
delimiter movement and small displacement is: 


y R.corr. rCorr. yn.apros. ; 


-1/ 


L + Cjjt 


^Corr.y R. 


aprox. 


(37) 


y R.Corr. TCorr. 


= <Pcorr.^^^^L + Cjjtsen 2 ^ 


EA 


t 


(38) 


2 1 EA’ I V M ^/L+C„t| 

Then the equation that models the force exerted at the ends or 
anchorages during the loading stage or deceleration of the impact body is: 


T^=AE 


YCorr. 2 




tsen 




2 

■Jf- 

(fr‘ 

+1-1 

VL+c„r 


(39) 


This equation is valid from the first contact up to a maximum tension is 
reached or only for the loading stage: max. ■ 

Hence, a model for the axial stress (as an average value for all the 
wires) at the ends is: 


(Tl=E 


<!> 


00 


Corr. 




tsen 




2 

■jf' 


+ 1-1 

^^L + CJJt 


(40) 


for I 

The last two equations represent a resonant system and only model the 
loading stage. This is why time is restricted up to the maximum tension is 
reached or the mass M is completely stopped. Then, just an instant later, 
the unloading stage begins. Thus, to model this process other equations are 
needed. To take into account the unloading stage, the elastic modulus of a 
wire rope has to be analyzed again. This parameter changes with every cycle 
and level of load or during unloading, and in general, in this stage its value is 
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bigger than when loading. Due to this dependency and knowing that a way to 
have a predictable design is having constant properties, there is a reason 
why it is recommended to load a new wire rope ten times before measuring 
its modulus of elasticity. But all these properties are given for quasi-static 
applications with little ratios of deformation f:10“^tol0“^fnfn/mm/s . And this 
is not enough to even classify a movement at the first scale of the 
Characterization of Weihrauch ( e with an order of 10°mm/mm/s ) [43], [69]. 
But high ratios of deformation are present during the transverse impact of 
wire ropes as the ones in the references used to validate the models here. 
There, the ratios can be up to 'e:2mm/mm/s . And this, of course, affects 
how the elastic modulus change when a wire rope is loaded. Observing the 
references results (used for validation) it can be say that as with a simple 
carbon steel bar, a high ratio of deformation makes it more resistant and 
changes the elastic modulus to a higher level than when it is loaded with an 
almost static load. The same is supposed for a wire rope and in this case its 
structure would be more stable. Above, unloading of a wire rope is described 
using the curve in figure 16, but every wire rope and specially a new one 
presents a permanent deformation after being loaded and unloaded (and 
more if there is plastic deformation). This is why a better approximation is 
given in figure 18, where a wire rope is unloaded from a maximum 
deformation to a final £^>0 with a higher elastic modulus, as it is near 
to be constant and more stable. 



Figure 18: Stress (a)-strain (e) curves of a 
loading and unloading cycle of an idealized 
wire rope. Embx, £f, are the maximum and 
minimum deformation respectively. 
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Also here is supposed that the unloading stage begins after the impact 
body reaches a null velocity and all its momentum has been transferred to 
the wire rope which, at the same time, has a maximum tension or a 
maximum deceleration. Then after this instantaneous tension the inverse 
begins to happen, this is, all the stress accumulated by the stress waves 
reflection begins to decrease in the same way as it increased but with more 
speed because the elasticity modulus is bigger. With geometry, the latter is 
expressed in figure 19, where an space without gravity or resistance medium 
is represented along with two perfectly rigid strings with uniform and identical 
properties are under longitudinal impact by a punctual mass which has an 
initial velocity Voo and leads to cycles of loading and unloading between a 
maximum and minimum tension (initial). This way if string A is unloading then 
string B is loading and when one reaches its maximum tension the other 
reaches its minimum. And because their properties are uniform then its points 
have the same deformation under the same tension. 

Figure 19 can have this interpretation; first, unloading of string A and 
loading of string B is expressed in the same time interval re[0;r„] , where 
string A or B reaches its maximum tension at time (the same for every 
point ) when loading or its minimum after unloading, this is, in this case 
tc_Mdx is how longer the loading stage last, and also coincides with the time 
that the unloading stage endures: to.Max ■ Now, if this time interval is divided 
into equal parts Ar=r„/n,then Ar=t,.-r,._i ; l<i<n, Vi,nelsi-{0} , thus 
it is possible to write the unloading tension of string A as a function of the 
tension in string B while it is loading for each time r, and every of its points 

and Sqjj . This is can be expressed as follows: 
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Ta{^o)—TB{^ n) y 
^a{^ 2)~^y 

TA{ti) = Ts{tn-i); 


And because t,= t,_i+At; then t„=t„_,.+ iAt and as iAt=t,- then 
t^_.=t^-ti, thus tension of string A while unloading can be expressed with 
r^(t.) = rB(r„-t;) and because the process is continuous, taking the limit, 
results in: 


T^{t) = TB{t„-t)=TB{-{t-t„)) (42) 

This is, the loading and unloading process are symmetrical functions of 

time with the values of tension transferred a time thus both stages 
take place at the same time. This is also compatible with the symmetry 
principle of time of Emy Noether, which says (basically) “ what is before is 


Unloading of string A Loading of string B 



Figure 19: Illustration to demonstrate the symmetry with respect to time of a 
loading and an unloading process of two identical perfectly elastic strings and 
with the same maximum and minimum tension levels T. Both are longitudinally 
impacted by a punctual mass M, t, represents equally spaced time instants, 
tn, represents the period from an minimum to a maximum and back to a 
minimum, g, is gravity, x,y, are the special coordinates, S, are the curvilinear 
coordinates. 


47 







after” and this way conservation of energy is established which in the case of 
strings under loading or loading is conservation of elastic potential energy 
[70]. Then, loading tension as a function of time and its symmetrical process 
of unloading of string A and B can be expressed in figure 20 as follows. 



Figure 20: Representation in a graphic Tension Vs Time of a loading and 
unloading process of two strings and according with the wire rope 
behavior of figure 19. Tmax, is the maximum tension while it is loaded or 
the beginning of the unloading stage, Tmm, this is the maximum tension at 
the beginning of the loading stage or when the unloading finish. To, level 
of tension while unloading de descarga, Tc.sim., level of tensions while 
unloading as a symmetric function of the loading stage. 

Now, the unloading tension To{t) of any ideal elastic string and with 

constant properties can be expressed as a symmetrical function of a loading 
process of the same string: 

^d( 0~ ^C.Sim.(“(^“^C.Mdx)) (43) 

This equation is valid at the interval tc.Mdx^^^^c.Mdx+^o.Mdx being the 
instant fc.Mdx > where the loading stage ends and the string tension reaches 
its maximum. But here (process of figure (20)), as both strings are the same 
then tc.Mdx=^D,Mdx ■ The same figure shows how a symmetric function can be 
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translated to any instant, hence, both stages can happen simultaneously by 
using expression -{t-tc_Mdx) and the unloading stage can happen just after 
the other is loaded using -(t-2tc.Mdx) ■ 

And because strings A and B have no mass or bending resistance, 
then, its state, in any moment, is equivalent to that of a string with a “delimiter 
movement”, with a triangular shape, small displacements and tension equal 
to the tension at the ends of a transversely impacted ideal string like in figure 
14. And this way is possible to continue approximating the transverse impact 
over wire ropes. 

But, as the change of modulus of elasticity after the loading stage is 
discontinuous this does not imply necessarily that changes in tension are 
also discontinuous. The maximum tension reached after loading 
Tc.Mdx = MAX{Tc) with from equation (39) has to be equal to the 
maximum tension reached after a unloading process symmetric to one of 
loading in the same wire rope with the unloading modulus, this is: 

Tc. Max ( . Max ) — Td. Mdx ( ^D. Max ) (44 ) 

In a real application the change in modulus is progressive but here its is 
possible to say that the unloading and loading modulus of elasticity are 
related with Ej^=RygE where i? 78 >l is a constant. And to show that 
equation (44) is true, equation (39) is used for both stages: 



As Ejj=R^gE , equation (46) can be expressed as follows: 
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( 47 ) 



In equation (47), , is the displacement during the unloading 

expressed as a symmetric function of the loading of the same string with its 
unloading modulus of elasticity and with the model of small 

displacements of equation (38). And because , the maximum 

displacement reached during the loading stage, is a constant then 
(yR.corr.(c.Ma./L/ 2 j ^=^79 's 3 constant too, then! 

■yj Ryg+i— 1 = ^1 Ryg‘^ Ryg + Ryg— Ryg (48) 

If roots in this relation are expanded using the first terms of a McClaurin 

series then the maximum tension reached is the same with which the 
unloading stage begins: 


1 "^79 

^78'*' ^78" 


’ ^79~ ^79 


(49) 


Knowing this, it is possible to express the model to approximate the 
tension at the ends during the unloading stage as a symmetrical function of 
the loading of the same wire rope but with its unloading modulus E^,: 


T D—AE j) 


^Corr.~ 




t rtFAir sen. 


} M 


V L+Cjjt^ 


(f)' 


+ 1-1 


(50) 


^ » — — 11 ~ j 

Valid in the time interval 1 c.max<^<^d.max+^c.max ■ Being, : 


^DESC — { {^d.mAx'^^c.mAx))) ’ ^c.mAx^^^^d.mAx'^^i 


C.MAX 


(51) 
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In (50) the symmetrical functions expressed as functions of {-t) have 
to be transferred to tc.MAx because it the unloading stage begins at that 
instant, and this is why tc.mAx'^'^D.MAx is added to -r. 


Before writing a model for the movement of the impact body, the model 
of the displacement during the unloading stage yR_corr.D of fh® delimiter 
string with small displacements has to be established. According to equation 
(50) displacement during the unloading stage as a symmetrical function of a 
loading process with module is: 


^00 \ M - I I M ""^"1 (52) 

2 Vl+c„wI 

And as is bigger than the one during loading then the maximum or 
initial value calculated with (52) is smaller to the final or maximum 
displacement reached during the loading stage: yn.corr.s{ED.Ma.)<yR.corr.(c.Ma.) - 
this is compatible with the minor deformation in the unloading stage showed 
in figure 18. Although, equation (52) does not model the displacement of the 
delimiter string middle point, it only indicates how much stretch has the string 
during the unloading stage. Then, to find out the actual displacement (small 
displacements model) in this stage in a way that the total displacement is a 
continuous function, the changes of coordinate of the middle point caused by 
its velocity has to be added, this is, the derivative of equation (52), 
y \ has to be integrated and added: 




E^A 

M 


*-DESC 




DESC 


10 In practice, determination of these times can be done with easiness using some software. 
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y R.Corr.D y R.Corr.{C.Mdx)'^ 


J v^d,, 


( 53 ) 


Vjj, is the velocity during the unloading stage, it goes from zero to a 
maximum when the string tension is zero again. This velocity is: 


Vn = 




Tr. 


Corr.S(ED) 


Where: 


' R.Corr.S{ED)~ y R.CorrJED 


)( (^D.MAX max))) 


(54) 


(55) 


And the displacement (small displacements model) during the 
unloading stage is: 


J' j' j ^ yR.Corr.S{ED)^t^ yR.Corr.S(ED)i^) Tfi. 


Corr. S{ ED) 


i^C.Mdx) 


(56) 


If time tc.Mdx is replaced into equation (55) then the term after sign 


minus in equation (56) is equal to the maximum reached during a symmetric 
loading process with an elasticity modulus equal to , as follows: 

y R. Corr. S{ED){^C. Max) — y R. Corr.{ED){^D. Max) — y R.Corr .{ED).MAX (57) 

And replacing this result into equation (53) then the displacement during 
the unloading is: 


.Corr.D ^ R. Corr .{C. Max) -Tr .Corr.{ED).MAx'^ y R.Corr.(ED)\ D.mAx"^^ C. mAx))) (^8) 

Valid in the time interval: 

^C.MAX^^ ^^D.MAX'*'^C.MAX (59) 

This way the model of small displacements is continuous in every 
instant no matter if the elastic modulus changes discontinuously. 

In sum, equations (38) and (58) model the small displacements of the 
middle point if an string which movement delimits that of an ideal string 
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impacted on its middle point and has the same tension than the ends of this 
ideal string. 

4.2 . Mathematical model of the impact body. 

Now, it is possible to model the kinematics of a body idealized as a 
material point that impacts an ideal string. This begins with the delimiter 
movement of a string with large displacements represented with the green 
string with triangular shape in figure 14. But first, as the lineal case has the 
following relation: 


F=2 


AEy, 


con -z-rr » 


(60) 


L/2 ’ -■ L/2 V 1 .V 4 +/, 

This means the angle a is assumed small and this way the stretch in 
the middle point So(L/2,r) can be approximated with the value . But here 


a better approximation is VTV4+y^-L/2 . It is important to emphasize that 
the tension of the ideal string is not that of the green string because the 
difference in shape and stretch. Then, and alternative to find a mathematical 
model is to establish that the stretch of the string with large displacements 
can be approximated for the loading stage y,^ ^ > ^i^cl unloading stage y,^ ^ , 
as a function of the displacements of the string with small displacements 
yn.corr. ■ Thus, for both stages: 


VLV4+/^-L/2«^VTV4 + ylc„,,-L/2 (61) 

With this, acceleration of point So{L/2,t) or of the impact body with 

mass M while the string is loading is: 




d,r 


c _ 2 

~ M 


Tr. 


Corr. 


Tr. 


Corr . 


I -/2 


Corr. j 


(62) 


for 0<r<rc.MAx 

And the acceleration while unloading is: 
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^ty M.D _ y R.Corr.D _ y R.Corr .D 

d ~ M L/2 J7^7d+v^ 

1 I ^ + yR.Corr.D 

^c.mAx^^'^^d.mAx'^^c.mAx 

Velocity of the impact body M in the string loading is: 


^tyM.C_^r r ( yR.Corr. 

d,t J ^ Ml L/2 Jl 




for 


Velocity while unloading is: 


^ty M.D. 


L 

- J 


^ ae ys 

.Corr.D 3^ R.Corr.D I W #■ 


Para el intervalo tc_MAx — t^^D.MAx'^^c.MAx 
In the equation above, Vp , is the velocity at the of the loading stage 

and is not necessarily zero in this model, because is approximated, this 
is,when the velocity is effectively zero, the calculated force is not necessarily 
maximum at the anchorages. Vp ^ can be calculated with: 


. MA 

=v - r 

^C. ^00 I 


^ AE ^ . Corr. _ yR .Corr. 1 W #■ 


for 

Displacement of the impact body M in the string loading is: 


t t 

s j 


^ AE Ir 

. Corr. . Corr. _I ^ 1 1 cf t 

M LI2 ' ' 


for 

Displacement while the string unloads is: 
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y M.D~^ 


F.C 


I. 
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\ 

y R. Corr .D 3^ i?. Corr. D ] 

d,t 

J 

J M 

^^2 ^1 L^/^+yl.Corr.ol 

.MAX 

^c.mAx 

1 



for Ic.mAx<^^^d.mAx'^^c.mAx 



d,t 


( 68 ) 


And finally the force of reaction F between the wire rope and the 
impact body is: 


And while unloading is: 


F^=-2AE^ 


. 1 . Corr. 

y R. Corr. \ 

1 L/2 

Vt / 4 + ^8 Corr. 1 

for 0^tKtc_MAx 

is: 


y R. Corr. D 

y R. Corr.D 

L/2 

Vt /4+ yii,Corr.D 


(69) 


(70) 


Para el intervalo tc_MAx — ^^^D.MAx'^^c.MAx 
Equations from (62) to (70) can be integrated numerically with a 

program like GNU Octave or using QtiPlot. 


Nd as a commentary, even with the discontinuity presented by the 
change of the elasticity modulus, all constructed here can be used directly as 
frontier conditions in a numerical method to solve the partial equations (8) 
and (9) with the only restriction that the string have to be elastic. This method 
can be similar to that developed by N. Cristescu for a string with a free end 
but now developing this is not at the reach of the author of this work [37], 


4.3 . Mathematical model of the internal stresses in 
the impact zone of a transversely impacted wire 
rope. 

The internal stresses arise due to the contact between the wire rope 
and the impact body. There emerges a force of reaction F, which is 
transferred to all the wires causing not only a resultant longitudinal tension 
but shearing, bending and torsion stress and also contact stress between 
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them. This way the wire rope presents higher stresses around the contact 
zone compared with what the anchorages support. Thus, the reason of this 
modeling is to include these additional stresses in the wires and then 
combine them with the longitudinal stresses of the wire rope and end up with 
a criterion to determine a security factor for the wire rope analyzed. 

To model these stresses there are some options, as: 

1. It is possible to consider the complete wire rope geometry and that of 
the impact body: the curvature of the surfaces in contact, of the wires, 
the force between each wire and other interior components. But this 
leads to a very elaborated geometry and represents another research 
work. 

2. Another alternative is to use only the surfaces in contact but this way 
the wires and strands re-arrangement is not taken in count. 

Then, a more viable option is to work as G. Rebel and R. Verreet did to 
study the surface wear of wire ropes under pressure when are rolled over 
drums. They modeled the contact stresses using Hertz theory, taking the 
nominal diameter of the wire ropes, a Poisson coefficient equal to 0.3 (as for 
steel) and the wire rope elastic modulus (smaller than that of steel) [74]. This 
way it is possible to approximate the stresses along the z axis(figures 21 
and 22), and only these ones because are the maximums around the contact 
zone and represent the additional stresses that arise when the wire rope is 
impacted^^ Two situations are considered: A steel cylinder with diameter Dj 
in contact with the equivalent cylindrical solid of the wire rope with diameter 
D 2 , their axis are at an angle of 90°; and the second case is a rigid plane 
front body of length n in contact with the wire rope equivalent cylindrical 
solid with diameter D. These cases are in concordance to what were used 


11 All the stresses in the contact zone and around it can be determined using the procedure and solving 
the elliptical integrals presented in references [71], [72], [73], 
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in the experiments and simulations used here to validate the models. The 
equivalent diameter of a wire rope is found using its metallic area: 


D^=D=2^[Ahc (71) 

The case of two perpendiculars cylinders is done assuming that the 
contact area is circular, next, an equivalent radio for the two cylinders in 
contact is found as follows [71], [72], [73]: 

(72) 

Where and are calculated using the diameters of both cylinders 
in contact y D^ . 


H, + H, = ^[2/D,+2/D,] 

H2-H^=^[4/Dl+2/Dl-8/{D^D^)Y^^ 
With (74) is calculated the contact area’s radio: 


(73) 

(74) 




11/3 


(75) 


In (75), fj, is the shear modulus E/2(1 +v). And the stresses along 

axis z are: 




1+ 


a 


(76) 


CTxx— Cfyy — Po 


1 2 -l\ 

(l+v)[^arctan(a/z)-l]+-(l+^) | 


(77) 


Po, is the maximum pressure reached over the circular contact area 


and depends on the maximum reaction force F 


MAX 


Pq—SFmAx/^^^ 


(78) 
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The equivalent perpendicular cylinders are as shown in figure 21: 



y| 


Figure 21: Cylinders of diameters Di, and D 2 in contact with its 
axis at 90°. x,y,z, are the coordinated axis with origin at the 
center of the contact zone. 

For a cylinder of diameter D in contact with a plane rigid body the 
contact is a rectangular area formed over the cylinder with a width b : 

j 1/2 

6 =(^FD[(l-v")/E]) (79) 

jtn 

The stresses along z axis in this case are: 


= -p„((l+2z^/fc^)(l+z^/fc^) ^'^-2z/b] 

(80) 

ayy=-2pQ v{{l+z^/bY'^-z/b) 

(81) 

a^=-p,{lWlbT'^ 

(82) 


Here, , is also the maximum pressure reached over the contact 
surface and depends on and on the contact length: 


_ mAx 

jrbn 

Both solids are represented in figure 22: 


(83) 
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Figure 22: A cylinder with radio R and a rigid plane of length h in 
contact. X, y, z are the coordinated axis with origin in the middle of the 
contact zone. 

Now it is possible to establish a model for the stress state around the 
contact zone during the transverse impact. This is represented in figure (23), 
there is included only the compression stresses which arise by the contact 
along z axis, and are also principal stresses cr^ , calculated for the 

maximum reaction force there is also added the stresses 

caused by the tension along the wire rope or the y axis, here 
<7 l.maa:(so=OjO which Can be calculated or measured. 



Figure 23: Maximum principal stress state along axis z inside a 
transversely impacted wire rope. Fmax ,is the reaction force 
between the impact body and the wire rope surface, Tmax. is the 
maximum tension reached by the wire rope at its anchorages, 
Oxx ,is the principal stress parallel to x axis, Oyy is the principal 
parallel to y axis, Ozz, is the principal stress parallel to z axis, a , 
longitudinal stress due to Tmax. 
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The foregoing lets to write a failure criterion for the impacted wire rope, 
that could be established using the Von Mises criterion (by analogy between 
a ductile material and a wire rope which structure rearrange while being 
loaded) to relate the internal stresses by contact and the longitudinal 
stresses, and then to compare results with the material strength limit and to 
determine a security factor. It has to be noted that the longitudinal stress at 
the point of impact (or zone) is not directly determined here but as the 
maximum value is needed, it is approximated using the maximum reached at 
the ends. Thus, using the Von Mises criterion along the z axis with origin at 
the impact point, the stresses are: 

((cTxx” CTyj" CTL.MAJf) ^L.MAX~ ^^zz) '^{o'zz~ (^xx) ) (^4) 

And this can be compared with the maximum strength of the wire rope 

^n.MAX ■ 

5 . Validation of the mathematical models. 

The experimental data used to validate the model of longitudinal forces 
and the movement of the impact body were collected in the experiments of C. 
Stolle and the simulated data were obtained in simulations done by Tytko et 
al y P. Gospodarczyk [31], [29], [33]. 

5.1 . Validation of the model of longitudinal stresses 
and the model of the impact body’s movement 
using the experimental data of C. Stolle. 

In this part, models are applied using properties and conditions as used 
by C. Stolle on its experiments. Tension at the ends of the wire rope while 
this is loaded or unloaded are calculated using equations (39) and (50), 
acceleration of the impact body is calculated with (62) and (63), its velocity 
with (64), (65), (66), and the displacement with equations (67) and (68). 
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C. Stolle did 6 experiments with steel wire ropes under transverse 
impact on its middle point, but here only experiments n° 1, 2, and 4^^ are 
taken. Wire rope characteristics necessary for mathematical calculations are 
in table 3. 


Wire rope type 

Number of strands 

Total number of 
wires 

Nominal 

diameter 

3x7 

3 

21 

19 mm 

Wire rope 
metallic area 

Weight per unit 
length 

Steel density 

Breaking 

strength 

1.54508E-04 

1.22803E+00 Kg/m 

7.80000E+03 Kg/m" 

1.80100E+05N 

Elastic limit 

Elastic modulus in 
an static test 

Elastic modulus of a 
pre-tensioned wire 
rope 

C. Stolle 
reco mended 
elastic modulus 

1.10400E+05N 

7.99500E+10Pa 

1.16000E+11 Pa 

1.04000E+11 Pa 


Table 3: Physical characteristics of a wire rope type 3x7 used by C. Stolle in 
his experiments [31]. 


The steel wire rope anchorages were conical sockets filled with epoxy 
resin and joined to a steel base which let the rotation of the anchorage 
around a vertical axis. This base, in turn, had a reinforced concrete base as 


shown in figure 24. 



Figure 24: Photograph showing the 
socket used as anchorage for the 3x7 
wire rope used by C. Stolle in his 
experiments [31]. 



Figure 25: Cross section of the 3x7 
wire rope used by C. Stolle in his 
experiments. 


12 C. Stolle, mentions that in his experiments 3,5 and 6 occur a lost of data and technical problems 
and the data were not appropriate for his simulation and are not for this work too. 
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C. Stolle, also built a car, specially done for his experiments, in its 
frontal part included a cylindrical steel tube of 076mm and A36 (figure 51) in 
a way that it was the only part in contact with the wire rope during the 
movement and no other part of it interfere with the car (impact body) 
movement. A rail was used to guide the tires and the car straight and 
perpendicular to the wire rope while it was impulsed by a pickup up to the 
wanted velocity of impact. Figure 26 shows the car contacting the wire rope: 



Figure 26: Impact position of the car and the 3x7 wire rope 

on the transverse impacting experiments of C. Stolle [31]. In this 

section the developed models are used and its results are compared with the 

experimental data of C. Stolle, and as him also did a numerical simulation 

this results will be used to compare the models. In his experiments the wire 

rope length was larger than the distance between anchorages: 46.2 mm 

more in experiments 1 and 2 and 50.8 mm in experiment 4. The wire rope 

was brand new and its shape before the impact had waves due to its winding 

over a drum and gravity. C. Stolle did included the additional length and the 

non-straight shape in its simulation and also added factors to model the 

internal friction, bending resistance, torsion and plastic deformation. 

But, in the other hand the mathematical models here do not take into 
account those factors. Then, the parameters to consider into the models are 
only the properties given in table 3, the car velocity at the instant of impact 
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Voo, its mass , the distance between anchorages L as represented in 
figure 27, and tables 8, 9 y 7: 


Rigid anchorages 


= [772.0; 772; 769.3] [ Kg] 

[6.625; 6.75; 10.63] [m/s] 
Wire rope 3x7. 



Figure 27: Idealization of experiments 1,2 and 4 of C. Stolle. L, is the 
distance between anchorages or ideal wire rope length. Me, is the 
impact body mass, Voo, is the velocity of impact, x,y are the coordinated 
axis and O, is the middle point of the wire rope and coordinate origin. 


5.1.1 . Modeling experiment 1 of C. Stoiie. 

Value of entry variables and the parameters considered here are in 
tables 4,6 and 5. Results are presented in figures 28, 29, 30 and 31 and 
tables 7,8,9 and 10. 


Impact velocity 

Mass of the car 

Distance 

between 

anchorages 

Initial wire rope 
tension (New wire rope 
without pre-tension) 

6.62500E+00 m/s 

7.72000E+02 Kg 

8.61100E+00m 

O.OOOOOE+OON 


Table 4: Entry variables to model experiment 1 of C. Stolle [31]. 


Impact body 
(car) mass M^, 

Impact body 
velocity 

Ideal cable length L. 

Elastic modulus 

E 

7.72000E+02 Kg 

6.62500E+00 m/s 

8.61100E+00m 

1.68411 E+11 Pa 

Wire rope 
metallic area 

Wire rope steel 
density 

Sound velocity in 
steel C|, 

Wire rope mass 

^cable 

1.54508E-04 

7.80000E+03 Kg/m^ 

4.64662E+03 m/s 

1.05746E+01 Kg 

M=M +M 



^c.MAX 

7.82575E+02 Kg 

1.81660E-02m 

3.64693E+02 m^ Vs 

1.84057E-01 Seg 


^D.MAX 

^C.MAX'*'^D.MAX 

1.44601E-01 Seg 

3.28658E-01 Seg 


Table 5: Parameters to model the experiment 1 of C. Stolle for the wire rope 
unloading stage. 















0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Time (s) 

Figure 28: Modeled tension T (KN) at the anchorages of the transversely impacted wire rope used in the experiment 1 of C. 
Stolle and his measured and simulated results. Tensions indicated with arrows are maximums. 


Max. Tension T as per 
the Math. Models. 

Max. Tension T 

Max. Tension T j 

1 i j u r.* II measured in the Exp. 

simu ated by C. Sto e. ^ „ 

^ 1 of C. Sto e 

Max. Calculated 
deformation ratio. 

1.56620E+05N 

1.61020E+05N 

1.65920E+05N 

8.40279E-02 mm/mm/s 

1.84057E-01 S 

1.94263E-01 S 

1.93208E-01 S 

6.69355E-02S 


Table 7: Comparison between the maximums of the tension at the anchorages of 
the transversely impacted wire rope of the Exp. 1 of C. Stolle. 
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Displacement Ym (m) 


0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Figure 29: Modeled displacement of the impact body with the measured and simulated displacement in the Exp. 1 of C. Stolle. 


Max. Displ. As per the 

Max. Displ simulated by 

Max. Displ. Measured in 

Math. Model y^. 

C. Stolle. 

the Exp. 1 of C. Stolle. 

9.75036E-01 m 

9.75400E-01 m 

9.75200E-01 m 

2.02940E-01 Seg. 

1.96260E-01 Seg. 

2.01420E-01 Seg. 


Table 8: Comparison of maximums displacements of the impact body, 
calculated with the measured and simulated in the Exp. 1 of C. Stolle. 
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0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Max. Return velocity 
as per the Math. Model 

V’m- 

,, _ ^ Max. Return velocity 

Max. Return Velocity .. .. _ ^ 

1 * -1 u 11 measured in the Exp. 

simulated by C. Stolle. ^ 

-5.30896E+00 m/s 

-4.25400E+00 m/s 

-4.54000E+00 m/s 

3.28090E-01 Seg. 

3.50000E-01 Seg. 

3.00000E-01 Seg. 


Table 9: Maximums return velocities of the impact body, calculated 
with the measured and simulated in the Exp. 1 of C. Stolle. 
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0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Time (s) 

Figure 31: Modeled acceleration of the impact body with the measured and simulated accelerations in the Exp. 1 of C. Stolle. 


Max. Acel. As per the Max. Acel. Simulated by 
Math. Model y”^ (g’s, c. Stolle (g’s, g=9.81 ml 
g=9.81 m/s^). s^). 

Max. Acel. Measured in 
the Exp. 1 of C. Stolle 
(g’s, g=9.81 m/s^). 

9.13274E+00g 

1.02500E+01 g 

9.49000E+00g 

1.84095E-01 Seg. 

1.71030E-01 Seg. 

1.93400E-01 Seg. 


Table 10: Maximums accelerations of the impact body, calculated 
with the measured and simulated in the Exp. 1 of C. Stolle. 
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5.1.1.1 .Discussion of the modeiing of experiment 1 
of C. Stoiie. 

Before using the developed equations, all the parameters have to be 
determined. Values of the material point mass M , its impact velocity Voo 
are not a problem, but as the wire rope has slackness between the 
anchorages, this requires a decision to make about what is the value of L, 
and here is taken as the anchorages distance knowing its difference with the 
real length is less than 50 mm. Another decision to make, and more 
important, is what value can be chosen for the elastic modulus knowing it 
changes with time and load; about this, C. Stolle recommends a value of 104 
Gpa for a new wire rope and 116 Gpa if it has a pretension but only to 
calculate the sound velocity C„ , being this a mean value because a new 
3x7 wire rope has a modulus of 79.9 Gpa on its first load test. C. Stolle does 
not give precisely what values he used to simulate the loading or unloading 
stage. Thus, in this point another decision has to be made about what value 
to use in the developed models, and knowing the aforementioned, a value of 
104 Gpa is used as a mean for the loading stage because the experimented 
wire rope is brand new. 

For the unloading stage the next reasoning is followed: finishing the 
loading stage its structure is more or less under a pretension, more 
stabilized, therefore its modulus would be approximately 116 Gpa, as it is 
recommended by C. Stolle, but, with the high deformation rate influence and 
as the elastic modulus in the unloading stage has to be bigger than in the 
loading stage then its modulus will change from 79.9 Gpa to 116 Gpa at the 
end of the loading stage and will change in the same proportion after 
beginning the unloading stage, this is, from 116 Gpa to 168.46 Gpa as is in 
tables 5 and 13. 

About the results, in figure 28, the mathematical models predict a higher 
tension during the loading up to after the referential elastic limit at -0.12 S 
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(being 110 KN, for an static test). This happens because the chosen modulus 
is constant and bigger than a new wire rope has at the beginning and also 
influence the wire rope pull through the socket’s resin which also absorbed 
energy and caused that a minor tension is needed to stretch it at the 
beginning as the experimental data shows, followed with a fast rise of tension 
after 0.05 S. Then its structure stabilize more and its tension reaches and 
surpass the elastic limit while its elastic modulus continues to increase, but, 
then the wire rope begins to deform plastically and decrease until the 
unloading begins. And is in these two parts where the approximation is better 
because the wire rope structure is even more stable but differs again, as the 
predicted tension is higher after the elastic limit and continues this way up to 
its tension becomes zero. This happens because, while unloading, the elastic 
modulus is not constant and decreases as shown in figure 10. 

Then, figure 31 shows how the calculated acceleration presents a 
marked discontinuity when the elastic modulus changes abruptly from 104 
Gpa, while loading, to 168 Gpa, while unloading. Before this jump the 
analytical curve can be considered as a mean value of the acceleration 
during the loading stage. While for the unloading stage this curve presents a 
better approximation, as also happens with tension, and allows predicting the 
maximum acceleration reached. This is an important value to consider when 
designing vial barriers because it is a reference to know the forces exerted 
over the passengers in a car when they have an accident. 

The predicted velocity in figure 30, approximates with more precision, 
both, the loading and unloading stages and even the final or return velocity. 
This happens because, even if acceleration has a large discontinuity, the 
experimental acceleration is not smooth, has peaks for some hundredths of a 
second that after integration results in a smooth velocity curve. And as the 
calculated acceleration is as a mean value during the loading stage and 
approximates better while unloading, this results in velocity values very close 
to the real ones. The predicted zero velocity does not coincides with the 
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maximum force, as was pointed out before; this is a result of the difference 
between tension, in value and time, between the ends and at the impact 
point. The same explanation has the displacement prediction showed in 
figure 29, which is almost the same in every instant except at the end. 

In general the analytical predictions are not as much close to the 
experimental data as are the C. Stolle simulations in every curve. But some 
disadvantages to use his work are that he does not provide all the necessary 
data to use his simulation, and LS-DYNA, is expensive and not libre 
software. However, the developed models here can be used if the elastic 
modulus is close to a constant, or in general if the structure properties are 
uniform or constants, this is, if mathematics is to be used then it is possible to 
stabilize the wire rope before being installed. This could decrease its capacity 
to absorb energy but a protection barrier has to have energy dissipators 
instead of allowing the wire ropes to have a plastic deformation near its 
breaking point or instead of building an over-engineered structure 
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5.2 . Modeling experiment 2 of C. Stoiie. 

Entry variables and parameters necessary to model this experiment are 
in tables 11, 12 and 13. Calculation results are presented in figures 32, 33, 
34 and 35 and in tables 14, 15, 16 and 17. 


Impact velocity 

Mass of the car 

Distance 

between 

anchorages 

Initial wire rope 
tension (New wire rope 
without pre-tension) 

6.75000E+00 m/s 

7.72000E+02 Kq 

8.61100E+00m 

O.OOOOOE+OON 


Table 11: Entry variables to model experiment 2 of C. Stolle [31]. 


Impact body 
(car) mass 

7.72000E+02 Kq 

Impact body velocity 

Voo. 

6.75000E+00 m/s 

Ideal cable length 
L. 

8.61100E+00m 

Elastic modulus E 

1.04000E+11 Pa 

Wire rope metallic 
area 

1.54508E-04 

Wire rope steel 
density 

7.80000E+03 Kq/m^ 

Sound velocity in 
steel C|| 

3.65148E+03 m/s 

Wire rope mass 

cable 

1.05746E+01 Kq 

M=M + M „ 

c cable. 


u. 

^c.MAX 

7.82575E+02 Kq 

2.35529E-02 m 

2.86589E+02 

1.83612E-01 Seq 

Table 12: Parameters to model the experiment 2 of C. Stolle for the wire rope 
loading stage. 

Impact body 
(car) mass 

7.72000E+02 Kq 

Impact body velocity 

Voo. 

6.75000E+00 m/s 

Ideal cable length 
L. 

8.61100E+00 m 

Elastic modulus E 

1.68410E+11 Pa 

Wire rope metallic 
area 

1.54508E-04 

Wire rope steel 
density 

7.80000E+03 Kq/m^ 

Sound velocity in 
steel C|| 

4.64661 E+03 m/s 

Wire rope mass 

^cable 

1.05746E+01 Kq 

M= M + M 

c cable. 

u. 


^c.MAX 

7.82575E+02 Kq 

1.85088E-02 m 

3.64692E+02 

1.83612E-01 Seq 


^D.MAX 

t +t 

*^c.MAX *^D.MAX 

1.43980E-01 Seg 

3.27592E-01 Seg 


Table 13: Parameters to model the experiment 2 of C. Stolle for the wire rope 
loading stage. 
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Figure 32: Modeled tension T (KN) at the anchorages of the transversely impacted wire rope used in the experiment 2 of C. 
Stolle and his measured and simulated results. Tensions indicated with arrows are maximums. 


Max. Tension T as per 
the Math. Models. 

Max. Tension T 
simulated by C. Stolle. 

Max. Tension T 
measured in the Exp. 

2 of C. Stolle 

Max. Calculated 
deformation ratio. 

1.62565E+05N 

1.60800E+05N 

1.65500E+05N 

8.72157E-02 mm/mm/s 

1.83612E-01 Seg. 

1.92900E-01 Seg. 

1.95610E-01 Seg. 

6.69355E-02Seg. 


Table 14: Comparison between the maximums of the tension at the anchorages of the 
transversely impacted wire rope of the Exp. 2 of C. Stolle. 
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Displacement ym (m) 


0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 



Figure 33: Modeled displacement of the impact body with the measured and simulated displacement in the Exp. 2 of C. Stolle. 


Max. Displ. As per the 
Math. Model y^. 

Max. Displ simulated by 
C. Stolle. 

Max. Displ. Measured in 
the Exp. 2 of C. Stolle. 

9.79615E-01 m 

1.03170E+00m 

1.05100E+00m 

1.99984E-01 Seg. 

1.99600E-01 Seg. 

2.08100E-01 Seg. 


Table 15: Comparison of maximums displacements of the impact 
body, calculated with the measured and simulated in the Exp. 2 of 
C. Stolle. 
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Time (s) 

Figure 34: Modeled velocity of the impact body with the measured and simulated velocity in the Exp. 2 of C. Stolle. 


Max. Return velocity 
as per the Math. Model 

y’M- 

^ w 1 * Max. Return velocity 

Max. Return Velocity .. .. _ ^ 

1 * j 1 - f-* 11 measured in the Exp. 

simulated by C. Stolle. ^ ^ 

-5.81938E+00 m/s 

-4.26130E+00 m/s 

-3.52000E+00 m/s 

3.27590E-01 Seg. 

3.00000E-01 Seg. 

2.81300E-01 Seg. 


Table 16: Maximums return velocities of the impact body, calculated 
with the measured and simulated in the Exp. 2 of C. Stolle. 
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Time (s) 

Figure 35: Modeled acceleration of the impact body with the measured and simulated accelerations in the Exp. 2 of C. Stolle. 


Max. Acel. As per the Max. Acel. Simulated by 
Math. Model y”^ (g’s, c. Stolle (g’s, g=9.81 ml 
g=9.81 m/s^). s^). 

Max. Acel. Measured 
in the Exp. 2 of C. 
Stolle (g’s, g=9.81 
m/s^). 

9.70000E+00g 

9.06940E+00g 

9.65377E+00g 

1.94500E-01 Seg. 

1.88190E-01 Seg. 

1.84094E-01 Seg. 


Table 17: Maximums accelerations of the impact body, calculated 
with the measured and simulated in the Exp. 2 of C. Stolle. 
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5.2.1.1 .Discussion of the modeiing of experiment 2 
of C. Stoiie. 

Here are used the same elastic modulus as for experiment 1 because 
the same new wire rope was used. Only the impact velocity changes a little 
percentage, from 6.625 m/s in exp.1 to 6.75 m/s, but this is enough to cause 
a larger difference between the calculated and measured tensions even after 
the elastic limit of 110 KN is reached while loading. But a better 
approximation of the maximum tension is also reached: in experiment 1 the 
model predicts 156.62 KN compared with a measured tension of 165.92 KN, 
while in experiment 2 the calculated tension is 162.56 KN with a measured 
165.5 KN. Although the mathematical curve has to have the same behavior 
than before (in figure 28) now these curves are closer to each other, and in 
general, to the simulated values. An explanation for all these effects is the 
larger stretch of the wire rope caused by the little increase of velocity which 
causes a bit larger plastic deformation and larger pull through the sockets 
than in Exp. 1: 1.5 mm in exp. 1 and 2.4 mm in exp. 2. This also means that 
a larger amount of energy was dissipated resulting in a minor measured 
tension as was observed by C. Stolle too, but not considered by the model, 
[31]. Figure 36, shows the wire rope pull through the sockets^^. 

The difference in tensions also causes an over-prediction of the 
maximum acceleration along the unloading stage, as showed in figure 35. 
And due to the bigger difference in acceleration, velocities also differ being 
lesser than the experimented while loading and greater when unloading (in 
module), then the displacement, figure 33 is smaller, specially, in the 
unloading stage. Also, as with experiment 1, the calculated zero velocity 
happens almost at the same time as the experimented but does not coincide 
with the maximum tension. 


13 It is worth noting here that the elastic wave refraction in the anchorages, which prevents reflection 
to duplicate the wave amplitude, can be included in the developed equations following the 
geometrical method presented in 4.1. 
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In sum, it is clear that when a wire rope is tensed up to its plastic zone 
and with a more changed elastic modulus and additional dissipation in its 
anchorages, then the mathematical model could predict the maximums but 
over-predict, with some percentage, tension level in both stages of the 
movement. 



Figure 36: Final pull through of 2.4 mm of the 3x7 wire rope in one of the 
sockets filled with resin used as anchorages in the experiment 2 of C. 
Stolle [31]. 


Something to note here is how the numerical simulation of C. Stolle also 
shows a greater difference with the experimental data, compared with his 
simulation of experiment 1. He also explain this with the larger permanent 
stretch of the wire rope and energy dissipation on the anchorages [31]. 
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5.2.2 . Modeling experiment 4 of C. Stoiie. 

Entry variables and parameters to model the experiment 4 of C. Stolle 
are in tables 18 and 19. Calculation results are presented in figure 37, 38, 39 
and 40 and tables 20, 21,22 and 23. 


Impact velocity 

Mass of the car 

Distance Initial wire rope 

between tension (New wire rope 

anchorages without pre-tension) 

1.06300E+01 m/s 

7.69300E+02 Kg 

8.61100E+00 m 

O.OOOOOE+OON 

Table 18: Entry variables to model experiment 4 of C. Stolle [31]. 

Impact body 
(car) mass 

7.6930E+02 Kg 

Impact body 
velocity 

1.05700E+01 m/s 

Ideal cable length 
L. 

8.61100E+00m 

Elastic modulus 

E 

1.04000E+11 Pa 

Wire rope 
metallic area 

1.54508E-04 

Wire rope steel 
density 

7.80000E+03 Kg/m 

Sound velocity in 
steel C|| 

3 3.65148E+03m/s 

Wire rope mass 

M K, 

cable 

1.05746E+01 Kg 

M= M + M „ 

u. 

u. 

^c.MAX 

7.79875E+02 Kg 

3.68184E-02m 

2.87084E+02 

7.93582E-02Seg 


Table 19: Parameters to model the experiment 4 of C. Stolle for the wire 
rope loading stage. 
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Time (s) 

Figure 37: Modeled tension T (KN) at the anchorages of the transversely impacted wire rope used in the experiment 4 of C. 
Stolle and his measured and simulated results. Tensions indicated with arrows are maximums. 


Max. Tension T as per 
the Math. Models. 

Max. Tension T 
simulated by C. Stolle. 

Max. Tension T 
measured in the Exp. 

4 of C. Stolle 

Max. Calculated 
deformation ratio. 

1.81000E+05N 

1.69300E+05N 

1.71260E+05N 

2.11925E-01 mm/mm/s 

7.33976E-02Seg. 

1.04000E-01 Seg. 

8.62000E-02Seg. 

6.62788E-02 Seg. 


Table 20: Comparison between the maximums of the tension at the anchorages of the 
transversely impacted wire rope of the Exp. 4 of C. Stolle. 
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Displacement ym (m) 


0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 



Time (s) 

Figure 38: Modeled displacement of the impact body with the measured and simulated displacement in the Exp. 4 of C. Stolle. 


Max. Displ. As per the 
Math. Model y^. 

Max. Displ simulated by 
C. Stolle. 

Max. Displ. Measured in 
the Exp. 4 of C. Stolle. 

Measured displacement 
when two strands break 

7.50116E-01 m 

1.06052E+00m 

1.13820E+00m 

8.79760E-01 m 

7.33976E-02Seg. 

1.08110E-01 Seg. 

1.19770E-01 Seg. 

8.62000E-02Seg. 


Table 21: Comparison of maximums displacements of the impact body, calculated with the 
measured and simulated in the Exp. 4 of C. Stolle. 
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Velocity y'w (m/s) 


0 0.05 0.1 0.15 



Time (s) 

Figure 39: Modeled velocity of the impact body with the measured and simulated velocity in the Exp. 4 of C. Stolle. 


Max. Return velocity 
as per the Math. Model 

_yV_ 

9.02364E+00 m/s 
7.33976E-02Seg. 


Max. Return Velocity 
simulated by C. Stolle. 

7.17060E+00 m/s 
1.08800E-01 Seg. 


Max. Return velocity 
measured in the Exp. 
4 of C. Stolle. 


Measured velocity when two 
strands break 


7.41000E+00 m/s 8.60550E+00 m/s 

1.19690E-01 Seg. 8.62000E-02 Seg. 


Table 22: Maximums return velocities of the impact body, calculated with the measured and 
simulated in the Exp. 4 of C. Stolle. 
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0 


0.02 0.04 0.06 0.08 0.1 0.12 



Time (s) 

Figure 40: Modeled acceleration of the impact body with the measured and simulated accelerations in the Exp. 4 of C. Stolle. 


Max. Acel. As per the Max. Acel. Simulated by 
Math. Model y”^ (g’s, c. Stolle (g’s, g=9.81 ml 
g=9.81 m/s^). s^). 

Max. Acel. Measured 
in the Exp. 2 of C. 
Stolle (g’s, g=9.81 
m/s^). 

Measured acceleration when 
two strands break (g’s, 
g=9.81 m/s^). 

7.04283E+00g 

1.10490E+01 g 

1.05960E+01 g 

1.04106E+01 g 

7.33976E-02 Seg. 

8.79130E-02 Seg. 

8.66130E-02 Seg. 

8.62000E-02 Seg. 


Table 23: Maximums accelerations of the impact body, calculated with the measured and simulated in the 
Exp. 2 of C. Stolle. 
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5.2.2.1 .Discussion of the modeiing of experiment 4 
of C. Stoiie. 

Here, only a loading stage is present, up to two strands break, this is 
why is used an elasticity modulus of 104 Gpa. Mathematical predictions of 
tension, acceleration, velocity and displacement have similar characteristics 
than those made for experiments 1 and 2: modeled acceleration is greater 
which results in a predicted velocity lesser than the measured while loading, 
although the displacement is very close up to the time when it breaks. 
Calculated tension is also greater because the model does not include the 
permanent stretch and the pull through the sockets. This is also why the 
predicted time when it breaks is less than the measured or simulated, then 
the maximum acceleration, velocity or displacement can not be predicted 
accurately. But is the steeped loading of the wire rope up to two strands 
break what provokes that the mathematical model can only predicts the wire 
rope breakage. C. Stolle could not approximate the behavior as well, and the 
basic reason is the no consideration of the wire rope geometry simplified to 
an uni-dimensional string or a set of points like C. Stolle does. Figure 57 
shows how the experiment finished with only one strand of the 3x7 wire rope 
without breaking but stretched. 

However, a measure to take in count is the level of tension at which the 
two strands break, 171.26 KN, while the maximum tension it can withstand is 
181 KN. But a difference of ~10 KN or 5.5% is not big for an engineering 
application because the security factor used for wire ropes is always more 
than 2.5 [75]. In addition, this confirm that assuming the anchorage tension 
as equal or similar to that of the impact point is very close to what reality 
shows and confirm that the wire rope inertia is not an important factor to 
consider when it is impacted with the aforementioned conditions. Although it 
does indicates a decrement in the wire rope breaking strength at the impact 
point and this is the main reason why the internal stresses model was 
developed. 
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5.3 . Validation of the model of longitudinal stresses 
and the model of the impact body’s movement 
using the simulation data of Tytko et al and P. 
Gospodarczyk. 

Both researches used a 1x7 wire rope, its characteristics are in table 24 


and figure 41: 


Wire rope type 

Number of strands 

Total number of 
wires 

Nominal diameter 

1x7 

1 

7 

15.7 mm 

Wire rope metallic 
area 

1.50000E-04 

Weight per unit 
length 

1.17200E+00 Kq/m 

Steel density 

7.80000E+03 Kq/m" 

Elastic limit 

2.46000E+02 N 

Breaking strength 

2.79000E+05N 

Elastic modulus of 
a wire rope with 
no pre-tesion 

1.95000E+11 Pa 




Table 24: Physical characteristics of the 1x7 steel wire rope used by Tytko et 
al and P. Gospodarczyk in its simulations [33]. 



Figure 41 : Cross section of a 1x7 wire rope 
with 015.7 mm. 

Both works used the software LS-DYNA (also used by C. Stolle), the 
authors do not give any detail about their numerical models, but they used 
the following conditions for the wire rope: An completely straight wire rope 
(simplified), its length was exactly equal to the distance between anchorages 
which were perfectly rigid. 
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Simulation of Tytko et al is modeled first, he does not present results for 
the kinematics of the impact body. Data and parameters to take in count are 
the wire rope properties of table 24, the block impact velocity Vjo, its mass 
and the distance between anchorages given in tables 25 and 26 and 
showed in figure 42. 

1 

=4000 Kg V= 3.0 m/s 

Rigid anchorages Wire rope 1x7 



Figure 42: Idealization of the simulation of Tytko et al. L is the wire rope 
length equal to the anchorages distance. Me, is the impact body mass, 
Voo, is the velocity of impact, x, is the coordinated axis and O, is the 
origin of coordinates or middle point of the wire rope. 

Tytko et al, presented the next figure 43 showing how his modeled wire 

rope breaks under the block after being impacted. 


jdLdeck by LS-P rePost 


Conto 

beam integration pt#1 
min=-0.678102, at ©tern# 2952 
max=1.95835, at eiem# 2970 




Figure 43: Transverse impact over a wire rope and its breakage 
simulated by Tytko et al using LS-DYNA. 

Below, the developed mathematical model is applied for the wire rope 
tension at the anchorages (or longitudinal) and compared with the simulation 
results obtained by Tytko et al. The work of Gospodarezyk requires another 
analysis this is why is given in the next section. 
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5.3.1 . Modeling the simulation of Tytko et al. 

Entry variables and parameters are in tables 25 and 26. Calculation 
results are presented in figure 44 and tables 27 [29]. 


Impact velocity 

Block mass 

Distance between 
anchorages 

Wire rope initial 
tension 

3.00000E+00 m/s 

4.00000E+03 Kg 

I.OOOOOE+OOm 

O.OOOOOE+00 N 

Table 25: Entry variables to model the simulation of Tytko et al [29]. 

Impact body 
(block) mass 

Impact body 
velocity 

Ideal cable length 
L. 

Elastic modulus 

E 

4.00000E+03 Kg 

3.00000E+00 m/s 

I.OOOOOE+OOm 

1.95000E+11 Pa 

Wire rope 
metallic area 

1.50000E-04 

Wire rope steel 
density 

7.80000E+03 Kg/m^ 

Sound velocity in 
steel C|, 

5.00000E+03 m/s 

Wire rope mass 

M K, 

cable 

1.17000E+00 Kg 

M=M +M 




4.00117E+03 Kg 

1.75437E-02 m 

1.71001E+02 m0 5/s 



Table 26: Parameters to model the simulation of Tytko et al for the wire 
rope loading stage up to its breakage. 
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Time (s) 

Figure 44: Modeled tension T(KN) in the 1x7 wire rope with the simulated tension of Tytko et al. 246 KN is the 
standard elastic limit, 279 KN is the standard strength limit and 289.63 KN is the simulated breaking tension. 


Time up to the elastic . t th i f i t Standard breaking Breaking strength Max. Deformation ratio 

limit as per the Math. ° ic imi strength used in the simulated by Tytko et calculated with the Math. 

I simulated by Tytko et al. .j ■ 

Model. ^ ^ Math. Model al. Model. 


1.99049E-02Seg. 

1.98825E-02 Seq. 

2.79000E+02 N 

2.89630E+05 N 

8.74344E-01 mm/mm/s 

2.46000E+05 N 

2.46000E+05 N 

2.12365E-02Seg. 

3.80749E-02Seg. 

2.12153E-02Seg. 


Table 27: Times to reach the elastic limit with the calculated and simulated breaking strength for the simulation of Tyko et al. 
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5.3.1.1 .Discussion of the modeiing of Tytko et ai’s 
simuiation. 

In this case only calculation of tension is necessary and then its 
comparison with the simulated data of Tytko et al. The closer approximation 
is due to the more ideal conditions imposed in Tytko et al’s numerical model: 
the anchorages are truly rigid and the wire rope length coincides with the 
distance between them, also the elastic modulus is 195 Gpa, which is very 
close to that of a steel bar: 207 Gpa, this is, is close to a constant. Hence and 
for these reasons, the mathematical model predicts accurately the simulation 
until the elastic limit of reference, 246 KN, where the slope of the numerical 
curve begins to decrease due to the plastic behavior of the wire rope until its 
breaks at 289.63 KN. This tension is actually greater than its static strength 
which does not consider hardening due to high deformation ratios. In the 
other hand, in the analytical model can be included, in posterior 
developments, the elastic modulus of the plastic zone and approximate better 
the simulation. 

An important observation here is how, even if the impact body is not a 
material point or does not have a narrow cylindrical front, instead it has a 
plane front, the mathematics is able to approximate the numerical result 
mainly because the large differences between the mass of the block and the 
wire rope. 
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5.3.2 . Modeling the simulation of Gospodarczyk. 

Entry variables to model this case are in table 28. 


Impact velocity 

Parallelepiped 

mass 

Distance between 
anchorages 

Wire rope initial 
tension 

7.50000E+01 m/s 

1.80000E+04 Kg 

7.50000E+00m 

O.OOOOOE+OON 

Wire rope additional 
lenqth 

Mass-wire rope 
contact shape 

Initial wire rope 
shape 

Number of wire 
ropes in the base 


Plane front of a 



O.OOOOOE+OOm 

parallelepiped (See 
figure). 

Straight 

420 


Table 28: Entry variables to model the simulation of P. Gospodarczyk [33]. 

Figure 45 shows the initial position, an instant before the impact 
between the 18000 Kg block, front of 2.5m x 1.5m at 75 m/s and an artificial 
base (barrier) that includes a sheet metal (with a non specified small 
thickness, the author does not indicates how much) with 420 1x7 wire ropes 
in a criss-cross pattern at 90° arranged in six levels 6*(2*13+4*11). The 
parallelepiped contact the base at its center. 



Figure 45: Initial position of the parallelepiped and base made of 
wire ropes an instant before the impact simulated by P. 
Gospodarczyk. The base has diameter of 07.5m and is used in 
deep shaft mines. 1 is the parallelepided, 2 is the set of wire 
ropes and 3 is a sheet metal [33]. 

The idealization begins supposing that only the wire ropes directly 
under the parallelepiped suffer the impact and that the sheet metal strength 
is negligible. This idea is shown in figure 46: 
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Figure 46: First idealization of the simulation of an impacted 
artificial base used in deep shaft mines and simulated by 
Gospodarczyk. 


Then, to be able to use the models, a second idealization has to be 
done, this way: the mass of the block is divided equally between all the wire 
ropes directly under it: 18000 Kg/126 = 142. 86 Kg per wire rope. And now 
this is the mass of the material point to be used into the models impacting a 
1x7 wire rope at 75 m/s as illustrated in figure 47. 


Figure 47: Second idealization of the simulation of an impacted artificial 
base used in deep shaft mines and simulated by Gospodarczyk. L, is the 
distance between anchorages that coincides with the diameter of the 
base. Me, is the mass of the block, Voo, is the impact velocity of the block, 

X is the coordinated axis and O, is the origin of coordinates and the middle 
point of the wire rope. 
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Then, as in the rest of cases, equation (39) is used, to calculate the 
tension. But here equation (29) of the model of small displacements is used 
to calculate the velocity of the block during the impact. All the reasons about 
why the equation of large displacements is not used are given in the 
discussion of this case in 5.3.2.1. Entry variables and parameters are in 
tables 28 and 29, calculation results are in figures 48, 49 and table 30. 


Impact body 
mass Mj. 

Velocity of impact 

Voo. 

Wire rope 
idealized length L. 

Elastic modulus 

E 

1.42860E+02 Kg 

7.50000E+01 m/s 

7.50000E+00m 

1.95000E+11 Pa 

Wire rope 
metallic area 

Wire rope steel 
density 

Sound velocity C,, 

Wire rope mass 

M M 

cable 

1.50000E-04 

7.80000E+03 Kg/m^ 

5.00000E+03 m/s 

8.77500E+00 Kg 

Mass in 

movement + 

^cable. 

u. 

U2 

tc.MAx@278.6 KN 

1.51635E+02 Kg 

8.53823E-02 m 

8.78402E+02 m° ^/s 

7.33144E-03Seg 


to.MAx@278.75 KN 

7.33462E-03 Seg 


Table 29: Parameters to model the simulation of Gosodarczyk for the wire 
ropes in an artificial base of a deep shaft mine during the loading stage. 
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Figure 48: Calculated T(KN) of the transversely impacted wire ropes in an artificial base used in deep shaft mines, 
simulated by Gospodarczyk. Tension indicated with an arrow is the breaking strength. 
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0 0.002 0.004 0.006 0.008 



Figure 49: Calculated velocity with the final velocity, from the simulation, reached by the 18 TN parallelepiped after passing 
through tha artificial base simulated by Gospodarczyk. Velocity indicated with and arrow is the calculated final velocity taking in 
count an standard breaking strength of 279 KN for a 1x7 wire rope. 


Final velocity as per the Math model (Small Final velocity simulated by 
displ. Model)@ 279 KN. Gospodarczyk. 

51.29 m/s 48 m/s 

0.0073388 Seg. 

Table 30: Final velocity calculated and simulated of the 18 TN block after 
passing through the artificial base simulated by Gospodarczyk. 
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5.3.2.1 .Discussion of the modeiing of 
Gospodarczyk’s simuiation. 

As here is used the same 1x7 wire rope as in the case of the work of 
Tytko et al, then the same elastic modulus is chosen, 195 Gpa, and also the 
same breaking strength (static) of 279 KN. The sole objective in this case is 
to find the final velocity of the parallelepiped after going through the entiry 
base breaking the sheet metal and the wire ropes. First, the time when the 
wire ropes break is calculated and with this is calculated the block velocity 
reached just after it breaks the last wire rope of the base. But to calculate this 
velocity is not a simply use of equations. Here the impact velocity of 75 m/s 
and the parallelepiped dimensions compared to the wire ropes length have to 
be considered before. The impact velocity corresponds in the third scale of 
the Weihrauch characterization, this is, plastic deformation is going to be 
present and up to the wire ropes break as is observed in Tytko et al’s 
simulation and with a small increase of tension. Also, the block plane front 
covers a considerable part of the wire ropes length of 7.5m: the 1.5m side 
covers a 20% and the 2.5m side a 33.3%. Therefore during the impact the 
wire rope points directly under the block only have deformations in an 
horizontal direction and its slope is going to be always zero or with a direction 
parallel to the plane surface of the rigid parallelepiped. The difference is how 
an ideal elastic string deforms when is impacted by a block with a shape or 
by a material point, in this situation the local inertia of the wire rope and its 
change driven by the stress waves provokes that the slopes of every point 
around the impact zone to be greater than the ones around the anchorages 
as is observed in figures 11 and 14, also N. Cristescu shows this on his 
numerical solution of the impact over a string with one end free and C. Stolle 
on the photographs of his experiments [31], [37], [53]. But these differences 
are not present when the parallelepiped impact the base, there, the points in 
the impacted zone stretches but do not contribute to the vertical 
displacement of the wire rope. 
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In sum, this is the answer to why is used the model of small 
displacements (light blue delimiter string) to approximate the results of the 
simulation with a predicted final velocity of 51.29 m/s very close to a 
simulated result of 48 m/s. While if the large displacements equations (62) 
and (65) are used then the final velocity predicted is slightly different to the 
impact velocity. The error in the approximation with the small displacements 
model happens because it does not include the sheet metal, any plastic 
deformation or any steel hardening by the high deformation ratios. Figure 50 
shows the simulation result. 



Figure 50: Final result of the simulation made by Gospodarczyk 
with the block decelerated from 75 m/s to 48 m/s [33]. 
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5.4. Qualitative comparison between the internal 
stresses model application and experimental 
and simulated observations. 

In this part, the experimental observations are the ones of C. Stolle and 
from the simulation only of P. Gospodarczyk’s work [31], [33]. Simulation of 
Tytko et al is not considered because they do not give the impact block 
measures. 

Here the model of internal stresses by the contact during the 
transverse impact over a wire rope is used. All the variable and parameters 
are used from the experiments and simulations specified above. Tension in 
the impact zone is assumed, by the reasons presented before, equal to the 
tension in the anchorages of the wire rope, in this case only the maximum 
value is used jointly with the maximum contact force predicted or 
measured (maximum acceleration multiplied by mass). Then, in order to use 
the model, a circular shape is assigned to the cross section of the wire rope 
equivalent to its metallic area as a sum of the cross sectional areas of each 
of its constituent wires. A maximum calculated acceleration is found from the 
maximum prediction during the loading or unloading using equations (69) and 
(70), the equivalent diameter is found with equation (72), the internal stresses 
along the transverse z axis of the equivalent perpendicular cylinders is 
found with equations (76), (77) and (78). Internal stresses for a plane rigid 
body impacting with a wire rope are found with equations (80), (81), (82) and 
(83). Finally an equivalent stress which combines the internal stresses in the 
contact zone with the longitudinal stresses are calculated using the Von 
Misses criterion with equation (84). 
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5.4.1 . Application of the internal stresses model to 
the experiment 1 of C. Stolle. 

Conditions and parameters of C. Stolle experiments are in 5.1. The 
experimental car he used had a front with a cylinder of 076mm made of A36 
steel as figure 51 shows [31]. 



Figure 51: Cylindrical impact front used in the experimental car 
ofC. Stolle [31]. 


Then the maximum contact force and tension for each case have to be 
calculated using the indicated equations. The 3x7 wire rope characteristics 
are in table 3 and figure 25. The same elasticity modulus used before to 
model the loading stage is used. A breaking strength cr^i^MAx is calculated 
with the breaking tension indicated by C. Stolle, 181 KN divided by the wire 
rope metallic area. Entry variables from experiment 1 of C. Stolle are in table 
31, calculation results are in figure 52 and table 32. 
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Wire rope metallic area 

1.5451E-04 m2 

Wire rope equivalent diameter 

1.4026E-02 m 

Cylindrical impact body diameter 

7.6000E-02 m2 

Wire rope Poisson ratio 

3.0000E-01 

Impact body Poisson ratio 

3.0000E-01 

Wire rope elastic modulus 

1.0400E+11 Pa 

A36 steel elastic modulus 

2.0700E+llPa 

Mciximum tension measured: 

Expl. Max. 

6.7395E+09 N 

Mciximum tension calculated: T 

mod. Mat.Expl. Max 

6.6538E+09 N 


Table 31 : Entry variables to apply the model of internal stresses 
to the experiment 1 of C. Stolle. 


Wire rope tension: Measured (Med.) and 
calculated (Cal.) 

T^^ =165.2 

T =156.62 KN 

Equivalent diameter of the cylinders in contact 

1.6325E-02 

1.6325E-02m 

Radio of the contact area: a. 

2.2719E-03 

2.2430E-03m 

Contact pressure Po. 

6.7395E+09 

6.6538E+09 

Maximum contact force measured in Exp.l 

F . ^ ,=M*a„. ^ , 

7.2855E+04 

7.0113E+04N 

Maximum longitudial stress 

1.0692E+09 

1.0137E+09 Pa 

Superficial Von- Mises stress 

2.0979E+09 

2.0365E+09 Pa 

Maximum internal Von Mises stress 

4.7992E+09 

4.7082E+09 Pa 

Minimum internal Von Mises stress 

V.Min 

1.0869E+09 

1.0307E+09 Pa 

Breaking strength 

1.1715E+09 

1.1715E+09 Pa 

% of the wire rope with stresses greater than the 
breaking strength over the axis Z of its cross 

43.00 % 

35.00 % 

sectional area. 




Table 32: Calculation results of the internal stresses model using data of the 
experiment 1 of C. Stolle. 
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Deep from the origin of the Z axis (% of the wire rope diameter) 

Figure 52: Results of the model of the internal stresses (maximums) around the impacted zone of the wire rope used in the 
experiment 1 of C. Stolle. With two curves one using the maximum calculated tension Tmod. Mat.Expi. Max, and the other with the 

measured tension TExpi.Max and compared with the wire rope breaking strength Qp.Max. 
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5.4.2 . Application of the internal stresses model to 
the experiment 2 of C. Stolle. 

Entry variables from the experiment 2 of C. Stolle are in table 33, 
calculation results are in figure 53 and table 34. 


Wire rope metallic area 

1.5451E-04 m2 

Wire rope equivalent diameter 

1.4026E-02 m 

Cylindrical impact body diameter 

7.6000E-02 m2 

Wire rope Poisson ratio 

3.0000E-01 

Cylindrical impact body Poisson ratio 

3.0000E-01 

Wire rope elastic modulus 

1.0400E+11 Pa 

A36 steel elastic modulus 

2.0700E+llPa 

Measured maximum tension: T,, , 

1.6550E+05 N 

Calculated maximum tension: T ^ 

mod. Mat.Exp2. Max 

1.6256E+05 N 


Table 33: Entry variables to apply the model of internal stresses 
to the experiment 2 of C. Stolle. 


Wire rope tension: Measured (Med.) and 
calculated (Cal.) 

T,,„, =165.5 

T a,c. =162.56 KN 

Equivalent diameter of the cylinders in contact 

1.6325E-02 

1.6325E-02 m 

Radio of the contact area: a. 

2.2885E-03 

2.2849E-03 m 

Contact pressure Po. 

6.7888E+09 

6.7780E+09 Pa 

Mciximum measured contact force in Exp.2 

F , „ ,=M*a„. „ , 

max.Exp2 Max.Exp2 

7.4467E+04 

7.4113E+04N 

Mciximum longitudinal stress 

1.0711E+09 

1.0521E+09 Pa 

Superficial Von Mises stress 0's Mm 

2.1084E+09 

2.0907E+09 Pa 

Mciximum Von- Mises stress 

4.8347E+09 

4.8157E+09 Pa 

Mciximum Von- Mises stress 

V.Min 

1.0893E+09 

1.0702E+09Pa 

Wire rope breaking strength o^,^^ 

1.1715E+09 

1.1715E+09Pa 

% of the wire rope with stresses greater than the 
breaking strength over the axis Z of its cross 
sectional area. 

44.00 % 

40.60 % 


Table 34: Calculation results of the internal stresses model using data of the 
experiment 2 of C. Stolle. 
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0% 20% 40% 60% 80% 100% 

Deep from the origin of the Z axis (% of the wire rope diameter) 

Figure 53: Results of the model of the internal stresses (maximums) around the impacted zone of the wire rope used in the 
experiment 2 of C. Stolle. With two curves one using the maximum calculated tension Tmod. Mat.Expi. Max , and the other with the 

measured tension Texpi. Max and compared with the wire rope breaking strength OR.Max. 
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5.4.2.1 .Discussion of the appiication of the internai 
stresses modei to the experiments 1 and 2 of 
C. Stoiie. 


Both predictions are similar. Using the measured data, results says that 
a 43% in exp. 1 and a 44% in Exp. 2 of the wire rope cross section in the 
impacted zone would to be broken. Whereas by using the calculated data the 
percentages decrease to 35% and 40% respectively. Although, nothing of 
these predictions happened, as are observed in figures 54 and 55. The wire 
ropes ended curved around the impact zone and none strand (which are a 
33.33% of the wire rope each one) broke. Thus, the model is not able to 


predict the experimental results in these cases. 



Figure 55: Permanent 
deformation of the wire 
rope around the impact 
zone after the Exp. 1 of 
C. Stolle [31]. 



Figure 54: Permanent 
deformation of the wire 
rope around the impact 
zone after the Exp. 2 of 
C. Stolle [31]. 


This error happens because the model does not take in count the 
rearrangement of wires and strands when are crushed by a transverse force. 
Stresses do not distribute as in a continuous cross section, its shape 
changes continuously while being impacted. The 3x7 wire rope has an 
arrangement showed in figure 25, each one of the three strands is held by 
the others only in two lines along a spiral, with no restriction to its lateral 
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movement'"^. This implies that there is not a continuous path where stresses 
can be transmitted all over the cross section as could happen to a larger 
extent in a 18x7 wire rope, for example, or in any other wire rope with more 
strands without a synthetic or natural core. 


5.4.3 . Application of the internal stresses model to 
the experiment 4 of C. Stolle. 

Entry variables to model the internal stresses are in table 35, calculation 
results are presented in figure 56 and table 36. As in experiment 4 of C. 
Stolle the wire rope breaks then only the loading elastic modulus is used and 
the same breaking strength are used as before in the other cases. In table 
37, are the calculations using two alternative tensions: one equal to the 
maximum measured and the other slightly higher. And are used with two 
maximum accelerations calculated accordingly (because the breaking 
strength changes then the time of rupture also changes). 


Wire rope equivalent diameter 

1.4026E-02 

m2 

Cylindrical impact body diameter 

7.6000E-02 

m2 

Wire rope Poisson ratio 

3.0000E-01 


Cylindrical impact body Poisson ratio 

3.0000E-01 


Wire rope elastic modulus 

1.0400E+11 

Pa 

A36 steel elastic modulus 

2.0700E+11 

Pa 

Maximum measured tension in Exp4: T„ „ 

1.7126E+05 

N 

Maximum breaking tension measured statically 
by C. Stolle: T ^ 

1.8100E+05 

N 

M^imum alternative tension: T . 

max. AIM 

1.7490E+05 

N 

Maximum calculated tension: T . 

max. Alt-2 

1.7126E+05 

N 


Table 35: Entry variables to apply the model of internal stresses 
to the experiment 4 of C. Stolle. 


14 A transverse force over a wire rope can also happens in other cases, for example a flattening can be 
caused by bending the rope severely over the rim of a sheave or any sharp object [76]. 
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Mciximum wire rope tensions: Measured Exp.4 
(Med.) and calculated (mcix. Strength) (Cal.) 

TMed. =171.26 

Tcaic. =181 KN 

Equivalent diameter of the cylinders in contact 

1.6325E-02 

1.6325E-02 m 

Radio of the contact area: a. 

2.3542E-03 

2.0545E-03 m 

Contact pressure Po. 

6.9837E+09 

6.0947E+09 Pa 

Mciximum measured contact force Exp.4 

F =M*a 

max. Exp2 Max. Exp2 

1.7126E+05 

1.8100E+05N 

Mciximum longitudinal stress 

1.1084E+09 

1.1715E+09Pa 

Superficial Von- Mises stress 

2.1743E+09 

2.0703E+09Pa 

Mciximum Von- Mises stress 

4.9478E+09 

4.4811E+09Pa 

Minimum Von- Mises stress 

V.Min 

1.1282E+09 

1.1845E+09Pa 

Wire rope breaking strength 

1.1715E+09 

1.1715E+09Pa 

% of the wire rope with stresses greater than the 
breaking strength over the cixis Z of its cross 
sectional area. 

58.00 % 

100.00 % 


Table 36: Calculation results of the internal stresses model using data of the 
experiment 4 of C. Stolle. 


Wire rope maximum alternative tensions 

=174.9 T, 

nsx. Alt-2. 171.26 KN 

Equivalent diameter of the cylinders in contact 

1.6325E-02 

1.6325E-02m 

Radio fo the contact area: a. 

2.0199E-03 

1.9988E-03m 

Contact pressure Po. 

5.9921E+09 

5.9293E+09Pa 

Maximum calculated contact force Exp.4 

F =M*a 

max.CalcExp2 Max.Calc.Exp2 

5.1205E+04 

4.9612E+04N 

Maximum longitudinal stress 

1.1320E+09 

1.1084E+09Pa 

Superficial Von- Mises stress 

2.0185E+09 

1.9873E+09Pa 

Maximum Von- Mises stress a,,,,. 

4.3916E+09 

4.3369E+09Pa 

Miminum Von- Mises stress 

V.Mm 

1.1444E+09 

1.1204E+09Pa 

Wire rope breaking strength 

1.1715E+09 

1.1715E+09Pa 

% of the wire rope with stresses greater than 
the breaking strength over the axis Z of its 
cross sectional area. 

58.00 % 

45.00 % 


Table 37: Calculation results of the internal stresses model using data of the 
experiment 4 of C. Stolle using two alternative maximum tensions. 
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Deep from the origin of the Z axis (% of the wire rope diameter) 

Figure 56: Results of the model of the internal stresses (maximums) around the impacted zone of the wire rope used in the 
experiment 4 of C. Stolle. With two curves one using the maximum calculated tension Tmod. Mat.Expi. Max , and the other with the 

measured tension Texpi. Max and compared with the wire rope breaking strength OR.Max. 
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5.4.3.1 .Discussion of the appiication of the internai 
stresses modei to the experiment 4 of C. 
Stoiie. 

Figure 56 shows how, when the wire rope tension at its anchorages is 
171.26 KN (maximum measured tension) then the model predicts that a 58% 
of the cross section has to be broken which almost coincide with the 66.66% 
that correspond to two strands broken as happened in the experiment 4 of C. 
Stolle. In the other hand if the maximum breaking strength is used, 181 KN 
all the wire rope have to be broken an this coincides with the supposition of 
the analytical model. And if the alternative maximum tensions are used with 
their correspondent calculated maximum accelerations then by using a 
tension of 171.26 KN (equal to the measured) then only a 45% would be 
broken but if an tension of 174. 9 KN is used then a 58% broken section is 
predicted and coincides with what is predicted using the experimental data. 
And is possible to make the prediction to coincide exactly with the reality 
observed of two strands, or a 66.66% broken section, by using a tension of 
176.5 KN and this would be the new referential and predicted breaking 
strength of the impacted wire rope. 

However, as calculations of experiments 1 and 2 are far from reality 
then the internal stresses model is inconsistent. And an explanation to why 
there is a coincidence with the results of experiment 4 is: first the 
rearrangement of wires and strands is done during the first milliseconds, 
next, a reconfigured cross section bears the impact force with small changes 
but as the impact velocity is 10.57 m/s, which is up to 56.6 % greater than the 
impact velocities of experiments 1 and 2, and with the tension reaching a 
maximum of 165.5 KN, then the breakage is inevitable. The particularity here 
is that the breaking tension was 171.26 KN which is not the maximum static 
measured of 181 KN and is also slightly bigger than the maximum tension in 
the other mentioned experiments. But at the end is enough to coincide with 
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the prediction of the internal stresses model. Only because the breaking 
strength is close to the achieved with a much less velocity. 

Figure 57 shows the final condition of the wire rope after finishing the 
experiment 4 of C. Stolle. Only one strand of the 3x7 wire rope stand 
unbroken and locally bent. 



Figure 57: final condition of the 3x7 wire rope after being impacted 
in the experiment 4 of C. Stolle. [31]. 
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5.4.4 . Application of the internal stresses model to 
the Gospodarczyk’s simulation. 

Characteristics of the wire ropes and the simuiation conditions are in 
5.3.2. The impact soiid has a piane front of 1.5m by 2.5m and faiis over wire 
ropes of 7.5m iong. Here, as in the simuiation, the soiid is considered rigid 
and its mass is divide between the 126 wire ropes (see 5.3.2.1) [33]. To use 
the modei is oniy necessary the finai contact force and the reached tension. 
Entry variabies are given in tabie 38, caicuiation resuits are presented in 
figure 58 and tabie 39. 


Wire rope metallic area 

1.5000E-04 

m2 

Wire rope equivalent diameter 

1.3820E-02 

m 

First side of the rigid parrallelepiped h-1. 

2.5000E+00 

m 

Second side of the rigid parallelepiped h-2. 

1.5000E+00 

m 

Wire rope Poisson ratio 

3.0000E-01 


Wire rope elastic modulus 

1.9500E+11 

Pa 

Maximum calculated break tension: T ^ -1 

mod. Mat.Gosp. 

side of 2.5m 

2.7875E+05 

N 

Maximum calculated break tension: T ^ -2 

mod. Mat.Gosp. 

side of 1.5m 

2.7856E+05 

N 


Table 38: Data from Gospodarczyk’s simuiation needed to use 
the modei of internai stresses. 


Calculated breaking tension for: side of 1.5m and 2.5m. 

Tca,c.i.5. =278.56 

"^0310.2.5^1 278.75 KN 

Half width of the contact area b. 

2.1989E-04 

1.7033E-04 m 

Contact pressure Po 

5.4554E+08 

1.1739E+09Pa 

Calculated mciximum contact force 

F =M*a 

max.Gosp. Max.Gosp. 

4.7108E+05 

4.7112E+05N 

Maximum longitudinal stress 

1.8583E+09 

1.8571E+09Pa 

Superficial Von- Mises stress 

2.0765E+09 

2.3266E+09 Pa 

Maximum de Von- Mises stress a„ 

2.0765E+09 

2.3266E+09 Pa 

Minimum de Von- Mises stress a„ 

V.Min 

1.8600E+09 

1.8600E+09 Pa 

Maximum breaking strength 

1.8600E+09 

1.8600E+09 Pa 

% of the wire rope with stresses greater than the 
breaking strength over the axis Z of its cross 
sectional area. 

100.00 % 

100.00 % 


Table 39: Calculation results of the internal stresses model using data of 
Gospodarczyk’s simulation. 
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Figure 58: Results of the model of the internal stresses (maximums) in the impacted zone of the wire ropes of the artificial base 
for deep shaft mines simulated by Gospodarczyk using two maximum calculated tensions Tmod. Mat.Expi. Max, one per each side of 

the parallelepiped that impact and compared with the wire rope breaking strength OR Max.- 
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5.4.4.1 .Discussion of the appiication of the internai 
stresses modei to the simuiation of 
Gospodarczyk. 

Prediction using the models indicates that the wire rope capacity due to 
a concentration of stresses at the impact (or contact) zone is not 
considerable. Statically measured the standard breaking strength of the wire 
rope is 279 KN then for the wire ropes under the 1.5m side of the 
parallelepiped the calculated breaking stress is 278.75 KN and for the ones 
under the 2.5m side is 278.56 KN. Why this happens has its base in the wire 
construction and properties. The arrangement of six wires around a central 
one prevents, to a larger extent compared with type 3x7 (for example), its 
relative lateral displacement keeping each other in contact. This kind of 
construction leads also to a high elastic modulus of 195 Gpa which is very 
close to that of a solid steel bar and means it does not change very much 
while being loaded or unloaded even for a brand new wire rope. And because 
of this the models can predict, at least qualitatively, what happens around the 
impact zone of the impacted wire ropes in the artificial base simulated by 
Gospodarczyk. 
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6 . Conclusions and recommendations. 

Transverse impact over steel wire ropes was investigated by building 
new mathematical models using mechanic’s general laws and then validating 
them with previous experimental and numerical works. 

The general objective of determining stresses along a wire rope at its 
middle point was partially achieved: 

✓ Two main models has been established: one of the reaction force at the 
ends or anchorages of the wire rope and another of the kinematics of 
the impact body. Both can be used either the wire rope have small or 
large displacements. A third model was done to analyze the decrement 
in strength due to the contact stresses that arise in the impact zone of 
the wire rope but, as is, it does not predict observed reality, although, is 
possible to improve it by choosing more appropriately wire rope 
properties The established models can predict with more accuracy if 
the anchorages used are rigid, the wire rope elasticity modules are 
constants or if it is stabilized and is made only of steel and round wires. 

✓ After assuming the wire rope as a cylindrical bar in contact with a 
deformable cylinder and a rigid parallelepiped, the internal stresses 
model based on a Hertzian contact stress theory was used to find the 
stresses around the contact zone of the impacted wire rope. But after 
comparing its results with experimental and numerical works, the model 
showed to be neither inconclusive nor appropriate basically because 
the parameters used like the Poisson ratio, do not represent or take in 
count all the behavior of the wire rope parts as its relative 
displacements or friction while being under the impact force. 

✓ Results of modeling the behavior of a simulated base made of wire 
ropes impacted on its center by a parallelepiped indicate that the 
models can serve as a basis to develop mathematical models of wire 
rope nets like the ones used in rock fall or landslide barriers. But more 
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experiments are needed in order to validate and know in deep the 
involved physics. 


6.1 . Conclusions of the validation of the model of 
longitudinal stresses and the movement of the 
impact body. 

Taking in count the validation using the experiments of C. Stolle and the 
simulation of Tytko et al, then the mathematical models can predict the 
anchorages tension of the wire rope and the movement of the body that 
impacts transversely it on its middle point. This happens as long as the 
anchorages are almost rigid and the wire rope modulus of elasticity is kept 
close to a constant during the loading and unloading stages or if the wire 
rope is stabilized before its use. Also, because the wire used by C. Stolle and 
Tytko et al are of different types, 3x7 and 1x7 respectively, then, at least for 
wire ropes with no natural or synthetic parts, helically stranded round profile 
wires and with the foregoing conditions, the established models are valid. 

And also because experiment 4 of C. Stolle and the simulation of Tytko 
et al showed that the wire rope breaking point under impact is around the 
impact zone and with at a tension close to its static breaking strength, this 
confirm the assumption that the anchorages tension is similar to the level of 
tension along the wire rope (under the established conditions). 

And regarding the modeling of the simulation of Gospodarczyk, it is not 
possible to say that the models can model the impact over bases made of 
wire ropes. Although, it is clear that the models approximate the block final 
velocity because the particular conditions of the simulation and specifically 
the near constant properties of the simulated 1x7 wire rope. 
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6.2. Conclusions about the application of the 
internal stresses model. 

Since predictions of observations in experiments 1 and 2 of C. Stolle did 
not coincide with the predictions done with this model, then it is not suitable 
for these cases. But this does not mean that is not appropriate for others 
because, basically, influence of the selected parameters was not investigated 
and only were chosen following the referenced work. An indication of this is 
the prediction done using data from experiment 4 of C. Stolle and the 
simulation of Gospodarczyk, this showed how the calculations coincided if 
the wire rope cross section does not change significantly during the impact. 
Despite the above, it is clear that a full metallic steel wire rope, in general, 
should have a Poisson ratio greater than 0.3 and also a transverse modulus 
of elasticity to be included in the model, hence, an orthotropic approximation 
can be better suited. And more, these results say that a model based in a 
Hertzian contact model can only works for strands and wire ropes with more 
than 3 strands and without natural of synthetic parts. 

Hence, more experimental test are needed in order to develop or 
improve the developed model and be able to predict the additional stresses 
that arise in the wire rope due to the contact with the impact body. 

6.3 . Summary of results. 

The next has been done: 

✓ By idealizing the wire rope as a string and modeling its transverse 
impact it was possible to calculate the tension or reaction force at the 
ends of the string or anchorages of the wire rope, that resulted to be 
similar to the tension at the impact point (middle point). 

✓ Also, by idealizing the impact body as a material point and simplifying 
its interaction with the wire rope, it was possible to approximate its 
kinematics. 
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✓ Validation of the mathematical models was carried out using 
parameters and results of previous numerical and simulation works that 
did not considered analytical approaches. This way, was showed that 
the model of the longitudinal stress and of the kinematics of the impact 
body can be used to predict the behavior of a transversely impacted 
wire rope. 

✓ As part of the mathematical models construction a so called method of 
“delimiter^® movements” was devised in order to have a path to follow 
when analyzing the movement of a transversely impacted string and 
then translated into a geometrical reference that could be translated, 
again, into mathematics, as the mathematical models of tension at the 
anchorages of the wire rope, kinematics of the impact body or of the 
internal stresses around the contact zone of the bodies. 

✓ The aforementioned method is of importance as it can be used in 
conjunction with the method of characteristics, with this can be 
approximated the border conditions for a string transversely impacted 
by a material point in a plane and then the strain, tension and 
kinematics of every point of the wire rope can be found by solving the 
hyperbolic partial differential equations. 

✓ A similar method was developed to model the unloading of an ideal 
impacted string, with which was demonstrated that this process is 
symmetric with respect to time of a loading process and expressed as 
the symmetric function (of time) of this process respectively. 

In sum, an engineering tool has been obtained with which is possible to 
design protection barriers and provide better and more solutions to vial 
accidents, fall of objects in mines, landslides and rockfalls. 


15 Bordering could be a better translation to “delimitante” 
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6.4 . Suggestions for future works. 

✓ In any other development or improvement to be done to the developed 
models, it should have its base in the philosophy of change, otherwise 
it could be another algorithm of solution and would not provide an 
insight about the phenomenon of transverse impact over wire ropes. 

✓ In the same way that here was developed and used the method of 
“delimiter movements” for an ideal string transversely impacted at its 
middle point, it is possible to apply the same method to model oblique 
impact and points others than the center. But before the cases with 
constant velocity must be studied, specially from this reference [35]. 

✓ When designing a protection barrier the security factor is going to be 
defined by the energy absorbers capacity to dissipate the impact, and 
how high are the additional stresses by contact in the wire rope. 

✓ A way to consider non rigid anchorages like the energy absorbers at 
the end of the wire rope is making a new geometrical analysis as 
presented in figure 17 and include the effect of refraction or a reflection 
that does not duplicate tension in equation (22). 

✓ Another suggestion, if the developed models are going to be used for 
design, then the wire rope answer have to be kept elastic. It is true that 
a software like LS-DYNA or another can be used but there is the work 
to find out what numerical algorithm and parameters are appropriate to 
be used in this kind of phenomena as the published works do not 
provide any code or the totality of the parameters they used. And if 
experiments are going to be done then the standard ASTM A741 -11 
(2016) have to be considered, there are suggested 3 types of wire 
ropes, but are generic denominations of various available constructions 
and this makes difficult to use only experimentation to design structures 
under impact like protection barriers. 
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8 . Annexes. 

8.1 . Annex 1. 

The aim here is to model the impact between a material point and a 
string from the first contact instant and before any wave reflects at its rigid 
ends. For this purpose figure 59 is presented: 



idealized curved line and its approximation with a triangular shape 
during its movement and before any waves reflects at its rigid ends. M, 
body mass; Voo, its initial velocity: V, it’s actual velocity; O, origin of axis 
X, y; L, is the distance between anchorage; F, represents the reaction 
force; Cn, is the longitudinal stress wave speed, X*, is the transverse 
wave speed and a is the angle between axis x and the triangular profile 
sides. 

Here is represented an ideal string during its first instants of movement, 
just after it is impacted by a material point with mass M and initial velocity 
. The transverse wave has not reached the ends yet and neither the 
longitudinal wave. The curved line is the idealization while the triangular 
profile with a angles is its approximation. With this is possible to formulate 
mathematical relations as follows: 

First, the reaction force is related with the string tension with: 

F=2Tsma (85) 

But tension T can be expressed with the deformation of the triangular 

string: 
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X. 

And then related with acceleration: 


( 86 ) 


F=M 


d.t^ 


(87) 


The extreme point x* of the triangle can be found with the transverse 
wave velocity: 


_ 

X = 




^ ,» 


^ cable P 


( 88 ) 


In equation (88), B , is a constant and T is the string tension and in a 
non linear relation: 



And replacing (88) in (89) an approximation for x* is obtained: 

* t 1/4 .- 3/4 -1 ,, 1/4 .B /2 1/2 /Qn\ 

X p «2 A p E t y (yU) 

Then, using equations (86) and (87) results, again, in an Endem-Fowler 

equation which models the movement of the material point: 




\ 1/2 / \ 1/2 
^ -J\ 1- 

.B -^13 B| 

I I 


(91) 


d/ M 

Now, by considering Yu. A. Dem’yanov’s supposition as in 4.1, the 
material point displacement has to be (like equation (11)) [66]: 


y(So=0,0 = yo(0=-Bit";Bi>0 (92) 

If this relation is replaced in equation (91), a solution is found when 

n=2 y B,=(Bj2f-. 


y{So=0,t) = {Bj2Ye 


(93) 
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Up to here everything seems to follow a logic path even all the 
approximations and simplifications done, but in coefficient of equation 
(91) the material point mass is dividing, indicating that while it increases then 
the material point will have less and less acceleration and as a consequence 
the reaction force F will decrease. However this has not any physical sense 
because if the body, keeping constant the velocity of impact, increase its 
mass then its momentum will increase too and after it contact an string with 
much more less mass and without bending strength, its points will accelerate 
more and with a bigger reaction force. In few words if the body increase its 
mass, it will sweep faster the string when impacting. 

An alternative is to use directly the general law of movement of Euler or 
general law of lineal momentum for an unidimensional body in a 
bidimensional space, this is, a string. Then the movement equation is: 

X* 

M^+J pA,,„,d,s^=MV (94) 

0 

This equation says: the impact body has a new velocity lesser than its 
impact velocity and with a new momentum that added to the string 
momentum is equal to the initial momentum of the impact body. However, 
this has to represent a wire rope being impacted and its impacted point has 
to accelerate from zero, this is, the first term to the left of the equal sign has 
to be zero but for an unidimensional string without bending strength and with 
a mass much small than the mass of the body, its point in contact with the 
material point will have the same velocity when contact each other in 
concordance with the assumption made in this thesis where the contacted 
string point instantaneously achieves a finite velocity after being impacted. 
And here is the contradiction, equation (94) should model a movement where 
the string point has to have a finite velocity but in reality the wire rope cross 
section or axis accelerates from zero when is impacted. Hence the behavior 
of the wire rope under transverse impact from the instant of contact to the 
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time when the impact body moves conjointly with the wire rope section over it 
impacts can not be modeled using an unidimensional string model and leads 
to a contradiction. Thus, here is an explanation: 

Idealizing a wire rope as an string and a solid as a material point 
provokes that after its impact a discontinuity on the tangent of the string. 
Nevertheless the Euler’s general laws of movement can only be used if the 
body considered is continuous and of class C3. And as the idealization made 
here uses an unidimensional string and a material point, it was possible to 
use Newton’s third law because it is only applicable to material points and not 
to continuities. Although Newton never formulated but only enunciated his 
laws, it was Euler who expressed mathematically for deformable continuous 
bodies [77], [78]. Citing the reference: 

"Every body perseveres in its state of rest or of uniform motion straight 
ahead, unless it be compelled to change that state by forces impressed upon 
it." In the generality maintained in modem mechanics, this axiom is not 
always valid, for a body may be subject to internal or external constraints not 
expressed in terms of a system of forces. For example, a rigid body subject 
to no applied force spins about some axis through its center of mass; its 
parts, which also are bodies, move in such a way that their centers of mass 
describe circles about that axis. Newton himself did not specify any 
mathematical properties of bodies or forces, and so his intentions must be 
inferred by the reader, and in the course of time different readers have read 
different meanings into his words.” [78] 

In conclusions the foregoing development also means that 
unidimensional idealizations can not be used to model the first instants of the 
impact between bodies (a time that can be called impact period instead to the 
whole interaction) because its not possible to take in count the relative 
velocity between its points in contact and up to them could share a common 
velocity. Furthermore, contact happens between surfaces and not necessarily 
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between the bodies gravity centers. This is why the usual momentum model 
of two spheres impacting can not consider tribological phenomena in and 
around the contact zone. In addition, momentum law used for material points 
only approximates what happens before and after the impact happens and 
therefore is no way to know the strain inside the impacted bodies. 

8.2 . Annex 2. 

Figure 60 shows the results of other 3 numerical models or codes 
developed for many years to study wire ropes under impact (longitudinal and 
transverse) and that were tried by C. Stolle in its thesis to compare the 
effectiveness of his own numerical model and simulation. In the next figure is 
also included the mathematical model of the anchorages tension applied to 
the experiment 1 of C. Stolle with the same intention of comparison. 
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Figure 60: Comparison of the mathematical model of the anchorages tension T(KN) applied to the transversely impacted wire 
rope used in the experiment 1 of C. Stolle, jointly with the simulation of C. Stolle and other numerical codes that he used to 
compare his results. [31]. 
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